Local nuclear energy density functional at next-to-next-to-next-to-leading order 
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We construct nuclear energy density functionals in terms of derivatives of densities up to sixth, 
next-to-next-to-next-to- leading order (N^LO). A phenomenological functional built in this way con- 
forms to the ideas of the density matrix expansion and is rooted in the expansions characteristic to 
effective theories. It builds on the standard functionals related to the contact and Skyrme forces, 
which constitute the zero-order (LO) and second-order (NLO) expansions, respectively. At N'^LO, 
the full functional with density-independent coupling constants, and with the isospin degree of free- 
dom taken into account, contains 376 terms, while the functionals restricted by the Galilean and 
gauge symmetries contain 100 and 42 terms, respectively. For functionals additionally restricted by 
the spherical, space-inversion, and time-reversal symmetries, the corresponding numbers of terms 
are equal to 100, 60, and 22, respectively. 

PACS numbers: 21.60.Jz, 21.30.Fe, 71.15.Mb 



I. INTRODUCTION 

The Skyrme force was introduced into nuclear physics 
more than half a century ago [1, 2] and it is still a con- 
cept that is widely used in methods striving to determine 
properties of nuclei irrespective of their mass number and 
isospin. However, at present we understand this concept 
in a significantly different way than it was originally pro- 
posed. Indeed, instead of using this force as an effective 
interaction within the Hartree-Fock (HF) approximation, 
we rather focus on the underlying Skyrme energy density 
functional (EDF), without direct references to the effec- 
tive interaction or HF approximation. 

In electronic systems, the use of functionals of den- 
sity is motivated by formal results originating from 
the Hohenberg-Kohn [3] and Kohn-Sham theorems [4], 
whereby exact ground-state energies of many-fermion 
systems can be obtained by minimizing a certain exact 
functional of the one-body density. This led to numer- 
ous extensions and applications, now collectively known 
under the name of density functional theory (DFT) [5-7] . 

The fact that properties of electronic systems are gov- 
erned by the well-known Coulomb interaction allows for 
derivations of functionals from first principles, by which 
token this approach can proudly be called a theory. For 
nuclear systems, the luxury of knowing the exact interac- 
tion is not there, so the analogous approaches developed 
in this domain of physics carry the name of the EDF 
methods. 

In this article we construct a phenomenological nuclear 
EDF based on strategies that are proper to effective the- 
ories [8]. There, guiding principles [9] are based on: (i) 
appropriate choice of effective fields, (ii) building effec- 
tive Lagrangian or Hamiltonian densities restricted only 
by symmetry principles, (iii) employing ideas of power 
counting. In the low-energy nuclear structure, correct 
fields can probably be associated with nonlocal one-body 
nuclear densities. Then, functionals of densities acquire 
the meaning of effective Hamiltonian densities. Although 



a formal construction of power-counting schemes is not 
yet available, ideas based on the density matrix expan- 
sion (DME) [10-16] (see also a recent example of an ap- 
plication to electronic systems in Ref. [17]) can be used 
to propose expansions in terms of moments of effective 
nuclear interactions, or equivalently, in orders of deriva- 
tives acting on the one-body densities. This is precisely 
the strategy we are going to follow in the present study. 

Effective field theories (EFT) were recently extensively 
applied in analyzing properties of nuclear systems. Here 
we are not able to give an even shortest possible review of 
this rapidly developing area of physics, but let us mention 
two specific examples. 

Firstly, the nucleon-nucleon (NN) scattering properties 
were very successfully described by employing the effec- 
tive nucleon-pion Lagrangian at next-to-next-to-next-to- 
leading order (N"^LO), see Ref. [18] and references cited 
therein. This showed that the EFT expansion is capa- 
ble of grasping the main features of nuclear interactions 
at low energies, without explicitly invoking microscopic 
foundations in terms of, e.g., heavy meson exchanges. 

Secondly, methods using the harmonic-oscillator effec- 
tive operators have been developed up to N'^LO, to be 
employed within the shell-model approaches [19]. There, 
the N'^LO expansion was explicitly expressed in the form 
of pseudopotentials that contain derivatives up to sixth 
order. Evidently, such pseudopotentials are exact equiv- 
alents of higher- order Skyrme-like forces. When averaged 
within the HF approximation, they would lead to EDFs 
depending on derivatives of densities up to sixth order. 
This allows us to label our approach with the traditional 
name of the N'^LO expansion too. 

There is also a recent significant effort in deriving nu- 
clear EDFs directly from low-energy QCD within chi- 
ral perturbation theory, see, e.g., Refs. [20-22]. This 
may have a potential of providing new important insight 
into the precise structure of terms in the EDF, while 
at present we are bound to proceed phenomenologically, 
with only the symmetry constraints available, as is done 



in the present study. 

In all rigorous EFT expansions, one strives to achieve 
convergence in describing physical observables by going 
to higher and higher orders of expansion. This is best 
illustrated by the so-called Lepage plots [19, 23], where 
for theories cut at different orders, relative errors of ob- 
servables are plotted as functions of energy. In nuclear 
EDF methods, this kind of convergence tests were never 
performed — simply because the functionals beyond the 
second order (NLO) of the standard Skyrme type had 
never been constructed or studied. The present study 
constitutes the first step towards this goal. 

Our paper is organized as follows. In Sec. II we de- 
fine basic building blocks for our construction, and then 
we construct local densities up to N'^LO. In Sec. Ill we 
construct terms in the EDF up to N'^LO and evaluate 
constraints imposed by the Galilean and gauge symme- 
tries. In Sec. IV we derive results for the case of con- 
served spherical, space-inversion, and time-reversal sym- 
metries. After formulating conclusions of the present 
study in Sec. V, in Appendix A we discuss general sym- 
metry properties of the energy density, in Appendix B we 
present details of the adopted choice of the phase conven- 
tion, and in Appendix C we list results pertaining to the 
Galilean and gauge symmetries. 



II. CONSTRUCTION OF LOCAL DENSITIES 

A. Building blocks 

Let p (rcr, r'cr') denote the one-body density matrix in 
space-spin coordinates. In what follows, in order to sim- 
plify the notation, wc omit the isospin degree of freedom, 
because in the particle-hole channel all densities appear 
in the isoscalar and isovector forms [24], and general- 
ization to proton-neutron systems does not present any 
problem. Within this assumption, the EDF we consider 
has the form: 



£ = / d?rnE{r), 



(1) 



where the energy density Ti.E{r) can be represented as a 
sum of the kinetic and potential energies. 



nEir)^—To+n{r). 
Zm 



(2) 



In the present study, we focus on the potential energy 
density Ti.{r) only. 

First, using the Pauli matrices CTq, where index a = 
{x, y, z} enumerates the Cartesian components of a vec- 
tor, the density matrix is separated into the standard 
scalar and vector parts [25], 

p {ra, r'a') = \p (r, r') S^a' + \'^{<jWa\ <j') Sa (r, r') , 



where 



(3) 
Pi^.r') = ^ p (rcr, r'cr) , (4) 

(J 

s{r,r') = ^p(ra,rV)(a'|^|a). (5) 



These two nonlocal densities will be used as building 
blocks of the functional together with the derivative op- 
erator V and the relative momentum operator fc. 



.^(V-V'). 



(6) 



To most easily satisfy the constraints imposed by the 
rotational invariance, in our method, the building blocks 
are represented as spherical tensor operators [26], i.e., 
P\^L {r, r') for A = and sa^ {r, r'), Va^, and /ca^ for A = 
1. In this notation, A is the rank of the tensor, and p = 
— A, . . . , +A is its tensor component. In the present study 
we use the following definitions of the building blocks in 
the spherical representation: 



Poo(r-,r') 
Si,f.={-i,o.i}{r,r') 

'^l,^l={-lfi,l} 



p(r,r'), 



{s^{r, r') - isy{r, r')) , s^{r, r'), -^ {s^{r, r') -I- isy{r, r'))^ , 



■ {;7| (V, - iVy) , V„ ^ (V, + iVy)} 



(7) 
(8) 

(9) 

(10) 



r 



In what follows, we most often omit indices and argu- In principle, arbitrary phase factors could be used in 

ments of these spherical tensors and we simply write p, front of the spherical tensors. In Appendix B, we discuss 
s, V, and k to lighten the notation. possible choices of such phase conventions, and determine 



the particular ones selected in Eqs. (7)-(10). These phase 
conventions, which are not the standard ones, are used 
throughout the paper and define the phase properties of 
all other objects that we construct by using the building 
blocks above. 



B. Higher-order derivative operators 

We begin by constructing all possible higher-order and 
higher-rank tensor operators from powers of the deriva- 
tive Vi^, where fi = — 1,0,+1 are the spin-projection 
components of the vector (rank-1) operator V. It is ob- 
vious that all possible nth-order powers of the derivative 
can be written as sums of terms Vi^^ . . . Vi^^. There- 
fore, any (n -I- l)th-order power is simply obtained by 
multiplying some nth-order power by a sum of Vi^ op- 
erators. Then, powers of a given rank can be obtained 
iteratively by vector coupling. 

In the second order, the two nabla operators can be 
coupled to angular momenta and 2. The coupling to an- 
gular momentum 0, [VV]o — A/VS, corresponds to the 
Laplacian operator. Furthermore, the coupling to angu- 
lar momentum 2, [VV]2, gives the second-order, rank-2 
derivative operator. The rank-1 coupling, [VV]i — 0, 
vanishes because the derivatives commute. Similarly, in 
each one higher order, a rank-L symmetric operator can 
be coupled with V only to L— 1 and L+ 1. Hence, all the 
nth-order powers have the form of A*-""^''^ multiplied 
by the _Lth-order rank-L (stretched) coupled operators 
for L == n, n — 2, . . . , (1)0. Then, up to N^LO, one ob- 
tains 16 different operators DnL listed in Table I. Any 
arbitrary tensor formed by coupled operators V can al- 
ways be rewritten as a sum of operators DnL through the 
repeated use of the 6j symbols. 

Exactly in the same way, we define 16 different oper- 
ators KnL, which are spherical tensors built of the rel- 
ative momentum operators k coupled up to N^LO, i.e., 
for n < 6 and L <Q. In the remainder of this section, we 
only discuss operators DnL , while all the results mutatis 
mutandis also pertain to operators KnL- 

The stretched coupled operators Dll for n = L, 



DLL = [y---[V[VV]2h...]L, 



(11) 



play a central role in our derivations below. They corre- 
spond to irreducible symmetric traceless Cartesian ten- 
sors built of the derivative V. They have 2L + I tensor 
components Dllm numbered by the quantum number 



J 



M = —L, . . . , L that we most often do not show below 
explicitly. Moreover, since terms in the EDF up to N'^LO 
depend only on operators Dll and Kll up to fourth 
order, L < A, see Sec. Ill A, below wc do not discuss 
stretched coupled operators of fifth or sixth orders. 

Equivalently, derivative operators Dll can be written 
in the Cartesian representation, in which their compo- 
nents are numbered by L Cartesian indices, DLL.ai...aLJ 
tti — x,y, z. The order of these indices does not matter 

TABLE I: Derivative operators DnL up to N'^LO as expressed 
through spherical tensor representation of the operator V. 



No. 


tensor DnL 


order n 


rank L 


1 


1 








2 


V 


1 


1 


3 


[vv]o 


2 





4 


[VV], 


2 


2 


5 


[VVloV 


3 


1 


6 


[V[VV],]3 


3 


3 


7 


[VV]^ 


4 





8 


[VV]„[VV1, 


4 


2 


9 


[ViVlVVlJg]^ 


4 


4 


10 


[VV]^V 


5 


1 


11 


[VV1JV[VV]2]3 


5 


3 


12 


[V[V[V[VV],]3]J3 


5 


5 


13 


[VV]^ 


6 





14 


[VV]^[VV]2 


6 


2 


15 


[VV]o[V[V[VV],]3], 


6 


4 


16 


[V[V[V[V[VVl,l3U5]g 


6 


6 



(totally symmetric tensors) and all traces vanish. 



E^^^>« 



0. 



(12) 



The Cartesian components DLL.ai...aL can be calcu- 
lated by using the detracer operator defined in Sec. 5 
of Ref. [27]. Up to fourth order they read: 

i^oo - 1, (13) 

Al,ai = V,,, (14) 

D22Mia2 — '^ai'^a2 ^ -gSaia2^j (15) 



J-^4:4,aia2a^a4 — Va^ V a2 ^^ 03 V 04 7 ^ i V a^ V a2 ^a3a4 1 VaiVa30a2a4 i vai^a4^a2a3 1 ^ a2^ as^aia^ 

'^ a2^ a^^aia^ ' ^a3 ^04^01 a2 ) ' 577^ I^aia2^a3a4 ^"*^ala3^a2a4 ' ^aia4^a2a^ 1 • \^ ' ) 



TABLE II: Spherical components of the derivative operators 
DnLM expressed through the Cartesian derivatives. Expres- 
sions for negative components can be obtained as DnL,-M = 
Klm, see Eqs. (B20) and (B22). 



-l)^-^'Dl, 



D„ 



Cartesian derivatives 



Duo = 
Dm = 

D200 = 

D220 = 

^221 = 

-D222 = 

-D330 = 

^331 = 

D332 = 
-D333 = 

-D440 = 






■^ 



^^ 



{dx + idy) dz 



»^^- 



^idl - -idl + 2dl) 



^(a 



,0 ,.. , idy) {dl + dl- Adl 

i\Vi{dx + idy) dz 



■■ {dx + id. 



1 



D441 = ^ (9:. + «9y) 9. (39^ + 39^ - 4af ) 



O442 = 
£'443 = 

D444 = 



% + 39,2 

-^{dx+idyf{dl^dl 



m 



We note here in passing that we could have equally 
well used the Cartesian derivative operators with traces 
not subtracted out, i.e., 



2?oo 


= 1, 


(18) 


^ll,ai 


- v„,, 


(19) 


^22,aia2 


- VaiVa2, 


(20) 


2^33, a 1020.3 


- Va,V„,Va3, 


(21) 


-^44,01020304 


- VaiVa2Va3Va4. 


(22) 



Representations (13)-(17) and (18)-(22) are equivalent 
in the sense that each operator -Dll. 01.. .01, is evidently 
a linear combination of operators A^^^^ ^^'^'DL>L'.ai...a^i 
for L' = L, L-2,...,(l)0. 

In principle, below one could replace the spherical rep- 
resentations of derivative operators shown in Table I by 
their Cartesian counterparts (13)-(17) or (18)-(22), and 
work entirely in the Cartesian representation. However, 
in our opinion, the use of the spherical representation is 
superior and more economical. Moreover, whenever cal- 
culation of the Cartesian derivatives is more suitable, we 
may express spherical components of the derivative op- 
erators through the Cartesian derivatives, as shown in 
Table II. An example of using the Cartesian representa- 
tion (18)-(22) is given in Sec. IV. 



TABLE III: Local primary densities (23) up to N-'^LO built 
from the scalar nonlocal density p{r,r') (v — 0). To simplify 
the notation the limit of r' = r is not shown explicitly. Stars 
(•) mark densities that enter the EDF up to N^LO. Bullets (•) 
mark densities that enter the EDF up to N^LO for conserved 
spherical, space-inversion, and time-reversal symmetries, see 
Sec. IV. The last two columns show the T and P parities 
defined in Eqs. (25) and (26), respectively. In addition, the 
time-even densities are marked by using the bold-face font. 



No. 


pnLvJ = 


density 


n 


L 


V 


J 


T 


P 


1 • • 


Poooo = 


[P]o 


















2 • 


piioi = 


[Ml 


1 


1 





1 


-1 


-1 


3 • • 


P2000 = 


[[fcfeloPlo 


2 















4 * • 


P2202 = 


[[kk],p]. 


2 


2 





2 






5 * 


P3101 = 


[[fcfeloMi 


3 


1 





1 


-1 


-1 


6 * 


P3303 = 


[[km^hph 


3 


3 





3 


-1 


-1 


7 * • 


P4000 = 


[[fcfcloP],, 


4 















8 * • 


P4202 = 


[[kk],[kk],p]. 


4 


2 





2 






9 


P4404 = 


llk[k[kk],]^],p]. 


4 


4 





4 






10 * 


psioi = 


[[fcfcJoMi 


5 


1 





1 


-1 


-1 


11 


p5303 = 


[[kkUk[kk]^]^p]^ 


5 


3 





3 


-1 


-1 


12 


P5505 = 


[[k[k[k[kk],]^],],p]^ 


5 


5 





5 


-1 


-1 


13 * • 


P6000 = 


[[fcfeloPlo 


6 















14 


P6202 = 


imim^p]^ 


6 


2 





2 






15 


P6404 = 


[[kk],[k[k[kk],],],p]. 


6 


4 





4 






16 


P6606 ~ 


iikik[k[k[kk],],],],],p]^ 


6 


6 





6 







C. Local densities 

Local densities are formed by acting several times on 
the scalar and vector nonlocal densities with the relative 
momentum operator k and taking the limit of r' = r. 
Using the spherical representation, the possible coupled 
fc-tensors (10) (up to sixth order in derivatives) Kn^ sue 
those given in Table I (replacing V with k). 

Acting with KnL on the scalar nonlocal density p (r, r') 
gives 16 different local densities up to N^LO (one for 
every term in Table I). They are listed in Table III. 
When acting with KnL on the vector nonlocal densities 
s (r, r'), one has to construct all possible ways of coupling 
the fc-tensors with the vector density. Obviously, each of 
the 4 scalar {L — 0) derivative operators gives one local 
density, while each of the 12 non-scalar (L > 0) derivative 
operators gives three local densities. Altogether, from the 
vector density one obtains 40 local densities up to N'^LO. 
They are listed in Table IV. 

Finally, all local densities can be denoted by four inte- 
gers nLvJ as 



PnLvj{r) = {[KnLPv{r,r')]j}^,^ 



(23) 



where the nth-order and rank-L relative derivative oper- 
ator KnL acts on the scalar {v — 0) or vector {v = 1) 
nonlocal density, and ranks L and v are then vector cou- 
pled to J. We call these local densities primary densities. 



TABLE IV; Same as in Table III but for densities built from 
the vector nonlocal density s (r, r') (v = 1). 



No. 


PnLvJ = 


density 


n 


L 




J 


T 


P 


17 * 


POOll = 


Ml 










1 


-1 




18 * 


Pino = 


[ks]. 


1 


1 









-1 


19 • • 


Pun = 


[ks]. 


1 


1 




1 




-1 


20 • 


P1112 ~ 


[ks]. 


1 


1 




2 




-1 


21 • 


P2011 ~ 


[[kk^s]^ 


2 







1 


-1 




22 • 


P2211 ~ 


mui 


2 


2 




1 


-1 




23 * 


P2212 ~ 


[[kk],s]. 


2 


2 




2 


-1 




24 * 


P2213 ~ 


[[kk]^s]^ 


2 


2 




3 


-1 




25 • 


P3110 ~ 


[[kk],ks]. 


3 


1 









-1 


26 • • 


P3111 = 


[[kk],ks]. 


3 


1 




1 




-1 


27 • 


P3112 ~ 


[[kk^ks]^ 


3 


1 




2 




-1 


28 • 


P3312 ~ 


[[k[kk],],s]. 


3 


3 




2 




-1 


29 • • 


P3313 = 


[[k[kk],],s]^ 


3 


3 




3 




-1 


30 • 


P3314 ~ 


mkkU,sl 


3 


3 




4 




-1 


31 * 


P4011 = 


imis]. 


4 







1 


-1 




32 • 


P4211 ~ 


llkkUkk],s], 


4 


2 




1 


-1 




33 • 


P4212 = 


[[kk],[kk],s]. 


4 


2 




2 


-1 




34 • 


P4213 ~ 


[[kkl[kk]^s].^ 


4 


2 




3 


-1 




35 • 


P4413 = 


[[k[k[kk],]^],s]^ 


4 


4 




3 


-1 




36 


P4414 = 


[[k[k[kk],],],s]. 


4 


4 




4 


-1 




37 


P4415 = 


[[k[k[kkU,ls], 


4 


4 




5 


-1 





TABLE IV: continued. 



No. 


pnLvJ = 


density 


n 


L 


V 


J 


T 


P 


38 * 


psiio = 


imiks]. 


5 


1 









-1 


39 • • 


P5111 = 


imiks]. 


5 


1 




1 




-1 


40 • 


P5112 = 


[[kkfoks]^ 


5 


1 




2 




"1 


41 • 


P5312 = 


[[kk],[k[kk],],s]. 


5 


3 




2 




-1 


42 


P5313 = 


[[kkl[k[kk]^y].^ 


5 


3 




3 




-1 


43 


P5314 = 


[[kkUk[kk],],s], 


5 


3 




4 




-1 


44 


P5514 = 


[[fe[fe[fc[fcfc]2]3]4]5*-]4 


5 


5 




4 




-1 


45 


P5515 = 


[[klk[k[kk]^]^iy]^^ 


5 


5 




5 




-1 


46 


P5516 = 


[[k[k[k[kk]^]^iy]^ 


5 


5 




6 




-1 


47 • 


P6011 = 


imu. 


6 







1 


-1 




48 • 


p6211 ~ 


imim.s]. 


6 


2 




1 


-1 




49 


p6212 = 


imim^s]^ 


6 


2 




2 


-1 




50 


P6213 ~ 


imim.s]. 


6 


2 




3 


-1 




51 


P6413 ~ 


[[kk],[k[k[kk],],],s]. 


6 


4 




3 


-1 




52 


P6414 = 


[[kk]^[k[k[kk]^]^ls]^ 


6 


4 




4 


-1 




53 


P6415 ~ 


[[kkUk[k[kk],]^],s], 


6 


4 




5 


-1 




54 


P6615 ~ 


[[k[k[k[k[kkU,],],],s] 


5 6 


6 




5 


-1 




55 


P6616 ~ 


[[k[k[k[k[kk]^]^l]^]^s] 


5 6 


6 




6 


-1 




56 


P6617 ~ 


[[klk[k[k[kk]^]^l]^]^s] 


.6 


6 




7 


-1 





The tensor components corresponding to the total rank 
J are not explicitly shown. 

One can also act on each of the local densities with 
derivative operators D„ii of Table I, and then couple 
ranks / and J to the total rank Q, i.e., 



PmI,nLvJ,Q{l') — [i'm/PriL^ j('")]q 



(24) 



For m > 0, we call these local densities secondary den- 
sities. We do not explicitly list them, because they can 
be obtained in a straightforward way from the primary 
densities corresponding to m = 0, PnLvj = Poo,inLvj,j, 
which are listed in Tables III and IV. 

In Tables III and IV, for completeness we also show the 
time-reversal (T) and space-inversion (P) parities defined 
as, 



T 
P 



(-1)". 



(25) 
(26) 



These definitions are based on the analysis of symmetry 
properties, which we present in Appendix A. To better 
visualize the time-reversal properties of the local den- 
sities, in Tables III and IV the time-even densities are 
shown in bold face. 

Local densities constructed above are complex. Tak- 
ing the complex conjugations gives relations derived in 
Appendix B: 



-M 



PmI,nLvJ,Q,M — \^^) Pml ,nLvJ,Q- M , 



(27) 



where the tensor components, denoted M, are shown 
explicitly. These relations allow for expressing positive 
tensor components through negative ones or vice versa. 
Therefore, complete information is contained in non- 
negative or non-positive tensor components only. The 
M — components are either real (for even Q) or imagi- 
nary (for odd Q), and hence 2Q + 1 real functions always 
suffice to describe a given local density of rank Q. More- 
over, all scalar densities are real, which was the basis of 
choosing this particular phase convention, as described 
in Appendix B. 



III. CONSTRUCTION OF THE ENERGY 
DENSITY 

A. Terms in the energy density 

Terms in the EDF we construct here are required to 
be quadratic in densities, invariant with respect to time 
reversal (Sec. Al), and covariant with respect to space 
inversion and rotations (Sec. A 2). All terms up to the 
N'^LO order in derivatives fulfilling these restrictions are 
constructed below. 

Using notation of Eq. (24), a general term in the energy 
density can be written in the following form. 



Tim I ,n L V J / \ r / \ 

mI,nLvJ,Q (^) = [Pm'/',n'L'«'J',Q (r )P: 



ml ,nLvJ,Q 



ir)]c 



(28) 



TABLE V; Numbers of terms defined in Eq. (28) of different 
orders in tiie EDF up to N'^LO. Numbers of terms depending 
on tlie time-even and time-odd densities are given separately. 
The last two columns give numbers of terms when the Galilean 
or gauge invariance is assumed, respectively, see Sec. IIIB. To 
take into account both isospin channels, the numbers of terms 
should be multiplied by a factor of two. 



order T-cvcn 


T 


-odd Total Galilean Gauge 


1 




1 


2 


2 2 


2 8 




10 


18 


12 12 


4 53 




61 


114 


45 29 


6 250 




274 


524 


129 54 



N-^LO 



312 346 658 



188 



97 



1000- 



u 

T3 



C 



T3 

O 



Eq. (30) 

Eq. (30), Galilean inv. 

Eq. (30), Gauge inv. 
•-•Eq. (28) 

▲-A Eq. (28), Galilean inv. 
■-■ Eq. (28), Gauge inv. 




100^ 



10^ 



2 4 6 

Order in derivatives 

FIG. 1: (Color online) Numbers of terms (28) and (30) shown 
in Tables V and VI, respectively, plotted in logarithmic scale 
as a function of the order in derivatives. 



where both densities must have the same rank Q to be 
coupled to a scalar. Moreover, their time-reversal and 
space-inversion parities (T and P) must be the same. 
Again, at N'^LO only terms with m' + n' + m + n < 6 are 
allowed. Then, the total energy density reads 



n{r) 



a 



'I ,n'L'v J rpra I' ,n' L' v J 
ml .nLvJjQ ml ,nLvJ.Q 



(r), (29) 



i.l.nLvJ,Q 



where C™j nLvj Q ^^^ coupling constants and the sum- 
mation runs over all allowed indices. 

Had we considered the case of coupling constants de- 
pending on density, all terms in Eq. (28) would have been 
independent of one another (up to a possible exchange of 
the two densities). Table V lists numbers of such inde- 
pendent terms, and they are also plotted in Fig. 1. 

In the present study, we concentrate on the case of 
density-independent coupling constants, in which case 
one can perform integrations by parts, so that the deriva- 
tive operators Dm'i' sue transferred from one density to 



TABLE VI: Same as in Table V, but for numbers of terms 
defined in Eq. (30). 

order T-even T-odd Total Galilean Gauge 






1 


1 


2 


2 


2 


2 


6 


6 


12 


7 


7 


4 


22 


23 


45 


15 


6 


6 


64 


65 


129 


26 


6 



N^LO 



93 



95 188 



50 



21 



the other. That this can always be done is obvious by 
the fact that the coupled derivative operators Dm'i' can 
always be expressed as sums of products of uncoupled 
derivatives Vi^ or V^. As a result of the integration by 
parts, the integral (28) can now be written as a sum of 
terms, where each term has the form: 



T 



l' L' v' j' 



mI,nLvj{r) = [Pn'L'v'.J'{r)[DmlPnLvj{r)].J']a, (30) 

where ranks / and J arc coupled to J'. Here, at N'^LO 
(i) only terms with n' + m + n < 6 are allowed, (ii) both 
densities must have the same time-reversal parity T, and 
(iii) their space-inversion parities must differ by factors 
(—1)^. Finally, in order to avoid double-counting one 
takes only terms with n' < n, and for n' — n only those 
with L' < L, and for L' — L only those with v' < v, and 
for v' = V only those with J' < J. Then, the total energy 
density reads 



n{r) = J2 



/-in'L,'V'J' rr\7l' L' v' J' f \ 

^"^ T '"T ■>, J^ rnl.nLvJV )i 



'i' L' v' j' 
'^ml ,nLvJ^ 



(31) 



where Cl^j'^j^'lj are coupling constants and the summa- 
tion again runs over all allowed indices. As we did for the 
local densities above, to better visualize the time-reversal 
characteristics of terms in the EDF, the coupling con- 
stants C^/'„i;^j corresponding to terms that depend on 
time-even densities are shown in bold face. 

Based on the results obtained in Sees. A 1 and A 2, and 
on Eqs. (25) and (26), we see that time- reversal invari- 
ance and space-inversion covariance require that 



_! -jji +v +n+v __ 1 



(32) 
(33) 



respectively. This means that integers v' -f w, n' -f n, 
and m must be simultaneously either even or odd. The 
numbers of all such allowed terms are given in Table VI 
and plotted in Fig. 1. The space- inversion covariance (33) 
requires that for all terms, the total orders in derivatives 
are even numbers, which defines our classification of the 
EDF up to LO (0), NLO (2), NNLO (4), and N^LO (6). 
In Appendix B, we presented terms in the EDF up to 
NLO, i.e., for zero and second orders, see Table XXll. 
The EDF at NLO is exactly equivalent to the standard 



TABLE VII: Terms in the EDF (30) that are of fourth or- 
der, depend on time-even densities, and are built from the 
scalar nonlocal density p{r,r'). Coupling constants corre- 
sponding to terms that depend on time-even densities are 
marked by using the bold-face font. Bullets (•) mark cou- 
pling constants corresponding to terms that do not vanish for 
conserved spherical, space-inversion, and time-reversal sym- 
metries, see Sec. IV. 



No. 


,— .n L V J 

^ml.nLv.l 


Pn'L'v'J' 


D^i 


pnLvJ 


1 • 


,-rOOOO 
*-^ 40, 0000 


[P]0 


[VV]^ 


Wo 


2 • 


<^20,2000 


[P]o 


[vv]o 


[[kk],p]. 


3 • 


,-.0000 
^-^22,2202 


[P]0 


[VV], 


[[kk]^p]2 


4 • 


^-,0000 
'-'00,4000 


Wo 


1 


[[kk]lp\^ 


5 • 


^2000 
'-'00,2000 


[[fcfeloPlo 


1 


[[kk],p]^ 


6 • 


,-»2202 
^-^00,2202 


\[kk]2p\2 


1 


[[kk]^P\2 



TABLE IX: Same as in Table VII but for terms that are built 
from the scalar nonlocal density p{r,r') and vector nonlocal 
density s{r, r'). 



No. 


^n'L'v'j' 
^mI,nLv,J 


Pn'L'v'J' 


DmI 


pnLvJ 


18 • 


^0000 
'-'31,1111 


Wo 


[VVloV 


[ksl 


19 • 


/-.OOOO 
'-'11,3111 


Wo 


V 


llkk],ks]i 


20 • 


/-,2000 
'-^11,1111 


[[fefeloPlo 


V 


[ks]. 


21 • 


y-,2202 
'-^11,1111 


[[kk]2p]2 


V 


[ks]. 


22 


/-,2202 
'-^11,1112 


[[kk\^p]2 


V 


[ksU 



TABLE X: Terms in the EDF (30) that are of fourth order, 
depend on time-odd densities, and are built from the scalar 
nonlocal density p{r,r'). 



TABLE VIII: Same as in Table VII but for terms that are 
built from the vector nonlocal density s{r,r'). 



Nr, rm'L'v'j' 
i\0. '~yrnI,nLvJ 



Pn'L'v'J' 



Dn 



PnLvJ 



No. 


r^n'L'v'j' 
'^ml.nLvJ 


Pn'L'v'J' 


Dn,! 


PnLvJ 


7 


,-,1110 
'-^20,1110 


[ks]o 


[VV]„ 


[ks]o 


8 


^1110 
'-'22,1112 


Mo 


[VV], 


[ks]^ 


9 


/-(lllO 
'-^00,3110 


[ks]o 


1 


[[fcfcloMo 


10 • 


^1111 
•^20,1111 


[ks]. 


[VV]„ 


[ks]. 


11 • 


^1111 
'^22,1111 


[ks]. 


[VV], 


[ks]. 


12 


/-(llll 
'-^22,1112 


[ks]. 


[VV], 


[ksU 


13 • 


,-,1111 
'-^00,3111 


[ks]. 


1 


[[fcfcloMi 


14 


^1112 
'-'20,1112 


[ks]^ 


[VV]„ 


[M2 


15 


/-(1112 
'-^22,1112 


[ksU 


[VV], 


[M2 


16 


,-,1112 
'-^00,3112 


[ks]^ 


1 


llkk]oks]2 


17 


/-<1112 
'-'00,3312 


[ks]^ 


1 


[[k[kk],^ls]^ 



Skyrme functional [28, 29], generalized to include all 
time-odd terms [24, 25, 30]. In both representations the 
functional depends, in general, on 14 coupling constants, 
and both sets are related by simple expressions given in 
Eqs. (B6)-(B19). 

In Tables VII-XVIII, we list all 45 and 129 terms in 
the EDF that are of fourth and sixth order, respectively. 
Together with 14 terms at NLO, listed in Table XXII, 
this constitutes the full list of 188 terms in the EDF at 
N^LO. 

After the complete list of terms in the EDF at N'^LO 
is constructed, one can check that not all of the local 
densities listed in Tables III and IV appear in the final 
EDF at N^LO. This is so, because it is not possible to 
couple all these densities to scalars, and simultaneously 
fulfill conditions (32) and (33), without obtaining more 
than total sixth order in derivatives. It turns out that 
out of the 56 local densities at N^LO, which are listed 
in Tables III and IV, only 35 occur in the final EDF 



23 


,-.1101 

C'20,1101 


[Ml 


[VV]„ 


[Mil 


24 


,-.1101 
L'22,1101 


[Ml 


[VV]^ 


[Mil 


25 


,^1101 
C'00,3101 


[Mil 


1 


[[fcfcloMli 



at N'^LO. In Tables III and IV such densities arc marked 
with stars (*) . Table XIX gives their numbers determined 
separately at each order. 



TABLE XI: Same as in Table X but for terms that are built 
from the vector nonlocal density s{r,r'). 



No. 


^.n'L'v'j' 
'-^ml.nLvJ 


Pn'L'v'J' 


D„,i 


PnLvJ 


26 


,-.0011 
^40,0011 


[sl 


[VV]^ 


[s]l 


27 


,-.0011 
L' 42, 0011 


[«]l 


[VVIJVV]^ 


K 


28 


,-.0011 
C'20,2011 


[«]l 


[VV]o 


im^s]. 


29 


,-.0011 
(-'22,2011 


[^]l 


[VV], 


mu. 


30 


,-.0011 
(-'20,2211 


[«]l 


[VV]„ 


im^s]. 


31 


,-.0011 
L'22,2211 


[s]. 


[VV], 


[[kkhs]. 


32 


,-.0011 
(-'22,2212 


[s]. 


[VV], 


[[kk]^s]^ 


33 


,-.0011 
C'22,2213 


[s]l 


[VV], 


[[kkhsh 


34 


,-.0011 
(-'00,4011 


[s]. 




imis]. 


35 


,-.0011 
(-'00,4211 


[s]. 




[[kkim^s]^ 


36 


,-.2011 
(-'00,2011 


l[kk]os]^ 




l[kk]os]i 


37 


,-,2011 
(-'00,2211 


[[kk]os]. 




[[kkhs]. 


38 


,-.2211 
(-'00,2211 


[[kkhs]. 




[[kk]^s]. 


39 


,-.2212 
(--00,2212 


[[kk]^s]^ 




[[kk]rs]2 


40 


,-.2213 
(--00,2213 


[[kkhsl 




[[kkhs]. 



TABLE XII: Same as in Table X but for terms that are built 
from the scalar nonlocal density p{r,r') and vector nonlocal 
density s{r,r'). 



TABLE XIV: Same as in Table XIII but for terms that are 
built from the vector nonlocal density s(r,r'). 



No. 


'^ml.nLvJ 


Pn'L'v'J' 


DmI 


pnLvJ 


41 


^1101 
<-^31,0011 


[Ml 


[VV]„V 


[s]. 


42 


^1101 
(-^11,2011 


[Ml 


V 


[[kk\Ai 


43 


^1101 

'-'11,2211 


[Mil 


V 


[[kk\Ai 


44 


^1101 

'-'11,2212 


[Mil 


V 


[[kk]A2 


45 


^3101 

'-'11,0011 


[[kk],kp], 


V 


\s], 



TABLE XIII: Terms in the EDF (30) that are of sixth or- 
der, depend on time-even densities, and are built from the 
scalar nonlocal density p{r,r'). Coupling constants corre- 
sponding to terms that depend on time-even densities are 
marked by using the bold-face font. Bullets (•) mark cou- 
pling constants corresponding to terms that do not vanish for 
conserved spherical, space-inversion, and time-reversal sym- 
metries, see Sec. IV. 



No. 


^mI,nLvJ 


Pn'L'v'J' 


D^i 


PnLvJ 


1 • 


,-,0000 
'-'60,0000 


[Plo 


[vvi^ 


[Plo 


2 • 


/-.OOOO 
'-'40,2000 


[Plo 


[VV]^ 


[[fcfcloPlo 


3 . 


^0000 
•-'42,2202 


[Plo 


[VVlo[VVl2 


[[kk],p]. 


4 • 


^0000 
'-'20,4000 


[Plo 


[vvio 


[[fcfcloPlo 


5 . 


/-.OOOO 
•-'22,4202 


[plo 


[vvi. 


[[fcfcljfcfc],pl. 


6 . 


,->0000 
•-'00,6000 


[plo 


1 


[[fcfcloPlo 


7 • 


^2000 
'-'20,2000 


l[kk]op]o 


[vvio 


l[kk]op]o 


8 . 


,->2000 
•-'22,2202 


[[kk]„p]. 


[vvi. 


[[fcfcl2pl2 


9 . 


/-.2000 
•-^00,4000 


mu]. 


1 


[[kk]lp]o 


10 • 


^2202 
'-'20,2202 


[[kk],p]. 


[vvio 


[[kk],p]. 


11 • 


,->2202 
•-'22,2202 


[[kk],p]. 


[vvi. 


[[fcfcl2pl2 


12 • 


/^2202 
'-'00,4202 


im^p]^ 


1 


[[kk],[kk],p]. 



B. Galilean and gauge invariance 

In the previous Sec. Ill A, the functional has been re- 
quired to be consistent with time reversal invariance, 
invariance under space reflections and rotational invari- 
ance. These constraints arise from symmetries of the 
NN interaction, see e.g. Refs. [31, 32]. Derivations above 
were much easier to perform in a general form, without 
imposing any other additional symmetry conditions. In 
this section, we treat such additional constraints coming 
from imposing the Galilean and gauge invariance. 

Assumption of the Galilean instead of Lorentz invari- 
ance goes hand in hand with using the Schrodinger equa- 
tion as a starting point and relies on the assumption 
that relativistic effects are negligible. This symmetry en- 
sures that the collective translational mass, calculated 
within the time-dependent HF or random-phase approx- 
imations, is correctly equal to the total mass, M — Am. 



No. 



/^n' L v' J' 
^ ml .riLvJ 



Pn'L'v'J' 



D„ 



PnLvJ 



13 

14 
15 
16 
17 
18 

19 • 

20 • 
21 

22 • 

23 • 
24 
25 

26 • 

27 • 



,^1110 
•-^40,1110 


[ks]. 


[vvi^ 


[Mo 


^1110 
•-^42,1112 


[Mo 


[vvioivvi^ 


[M2 


•-^20,3110 


[Mo 


[vvio 


[[MoMo 


,^1110 
•-^22,3112 


[Mo 


[vvi. 


[[M0M2 


^1110 
•-^22,3312 


[Mo 


[VV], 


[[k[kkU-A2 


^1110 
'-'00,5110 


[Mo 


1 


[[MoMo 


,^1111 
•-^40,1111 


[Ml 


[VV]^ 


[Ml 


^1111 
•-^42,1111 


[Ml 


[vvio[vvi2 


[Ml 


•-^42,1112 


[Ml 


[vvio[vvi2 


[M2 


'-'20,3111 


[Ml 


[vvio 


[[MoMi 


•-^22,3111 


[Ml 


[vvi. 


[[MoMi 


/^llll 
'-'22,3112 


[Ml 


[vvi. 


[[M0M2 


'-'22,3312 


[Ml 


[vvi. 


[[MM2I3SI2 


^1111 
•-^22,3313 


[Ml 


[vvi. 


[[k[kkUA, 


'-'00,5111 


[Ml 


1 


[[MnMi 



Therefore, in principle, the Galilean invariance should al- 
ways be imposed. However, in many phcnomenological 
approaches, like the non-interacting or interacting shell 
model, the Galilean symmetry is not considered, because 
the translational motion is not within the scope of such 
models. The question of whether the Galilean symmetry 
must be imposed in phcnomenological models is not yet 
resolved, and in the present study we keep this question 
open. 

For a local interaction, f (rj,r2,ri,r2) = 5{r'i — 
ri)5{r2 — T2)f(ri,r2), the HF interaction energy is in- 
variant with respect to the local gauge [36]. Therefore, 
for the total energy (1) obtained within the EDF method, 
one may also consider constraints resulting from assum- 
ing local gauge invariance. As mentioned, this symmetry 
is only fulfilled when the forces involved are local. An 
example of a local approximation is the well known lo- 
cal one pion-exchange (OPE) potential, which is only an 
approximate representation of the correct nonlocal OPE 
Feynman amplitude [33]. This approximation is good 
as long as the relative momenta of interacting particles 
are about the same in the initial and final states (see 
Fig. 10 in Ref. [33]). Some of the fitted NN potentials 
like the Argonne Vis, Nijm-II, and Reid93 use this local 
approximation while others (CD-Bonn) use the full non- 
local OPE amplitude [33]. The most important nonlocal 
term is the two-body spin orbit interaction [32], which 
violates the assumption of gauge invariance. However, 
in this case a gauge invariant spin-orbit term (used in 
the Skyrme and Gogny forces) can be obtained in the 
short-range hmit [32, 34, 36]. 

For the EDF derived in this work it is, however, the 
symmetries of effective forces rather than bare forces that 
should be considered. One of the methods to obtain an 



TABLE XIV: continued. 



No. 






n L V J 



Pn'L'v'J' 



Dr, 



PnLvJ 



28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 • 

43 

44 

45 

46 • 

47 



^1112 
<-^40,1112 


[ks]^ 


NWo 


[ks]2 


,^1112 
*-^42,1112 


[ks]. 


[VVIJVVI2 


[ks]^ 


<-^44,1112 


[ks]. 


[V[V[VV]J3]4 


[ks]^ 


<-^22,3110 


[ks]. 


[VV]^ 


[[kk],ks]. 


^1112 
'-^22,3111 


[ks]. 


[VV]^ 


{[kk]oks]. 


<-,1112 
'-^20,3112 


[ks]. 


[VV]o 


im.ks]. 


y^lll2 

'-'22,3112 


[ks]. 


[VV]^ 


[[kk],ks]. 


^1112 
'-'20,3312 


[ks]. 


[VV]o 


[[k[kk],]^s]. 


<-,1112 
'-'22,3312 


[ks]. 


[VV]^ 


[[k[kk],],s]. 


y^lll2 

'-'22,3313 


[ks]. 


[VV], 


l[k[kkU,s]^ 


^1112 
'-^22,3314 


[ksh 


[VV]^ 


mkkU.s], 


<-,1112 
'-^00,5112 


[ks]. 




imlks]. 


/^1112 
'-'00,5312 


[ks]. 




[[kkUk[kk]2]^s]2 


<-,3110 
'-^00,3110 


[[kk],ks]o 




[[kk]oks]. 


<-,3111 
<-^00,3111 


[[kk]oks]. 




imMi 


/^3112 
'-'00,3112 


[[kk],ks]. 




[lkk]M2 


<-,3112 
'-'00,3312 


[[kk],ks]2 




[lk[kk],],s]. 


,-,3312 
'-'00,3312 


[[k[kkU,s], 




[lk[kkU,s], 


y^3313 

'-'00,3313 


[[k[kk]^],s]. 




[[k[kkU,s], 


,-,3314 
'-^00,3314 


[[klkk],],sl 




l[k[kk]2]^sl 



TABLE XV: Same as in Table XIII but for terms that are built 
from the scalar nonlocal density p{r,r') and vector nonlocal 
density s{r,r'). 



No. 



> r / / T> 

y^n Li V J 

^ml ,nLvJ 



Pn'L'v'J' 



D^ 



pnLvJ 



48 • 


,-,0000 

'-'51,1111 


Mo 


[VV]^V 


[ksh 


49 • 


/-.OOOO 
'-'31,3111 


Wo 


[VV]oV 


im.ksi 


50 • 


,-,0000 
'-'33,3313 


Mo 


[V[vvi,i3 


mikkU^s]^ 


51 • 


^0000 
'-'11,5111 


[P]0 


V 


imlks]. 


52 • 


/-,2000 
'-'31,1111 


[[fcfcloPlo 


[VV]„V 


[ks]. 


53 • 


/^2000 
'-'11,3111 


[lkk]oP]o 


V 


[[kk]gks]^ 


54 • 


,-,2202 
'-'31,1111 


[[kk],ph 


[VVloV 


[ks]. 


55 • 


,-,2202 
'-'33,1111 


[[kk],p]. 


[V[vvi,i3 


[ks]. 


56 


/-,2202 
'-'31,1112 


[[kk],p]. 


[vvIqV 


[ks]. 


57 


,-,2202 
'-'33,1112 


[[kk],p]. 


[V[vvi2i3 


[ksh 


58 • 


,-,2202 
'-'11,3111 


{[kk],p]. 


V 


[[kk^ksl 


59 


,-,2202 
'-'11,3112 


[[kk],p]. 


V 


[[kk],ks]2 


60 


,-,2202 
'-'11,3312 


[[kk],ph 


V 


[[k[kkU,s], 


61 • 


,-,2202 
'-'11,3313 


{[kk],p]. 


V 


{[k[kk],],s].^ 


62 • 


,-,4000 
'-'11,1111 


imu. 


V 


[ks]. 


63 • 


,-,4202 
'-'11,1111 


[[kk\,[kk],p]. 


V 


[ksl 


64 


,-,4202 
•-'11,1112 


[[kkUkk],p], 


V 


[ks]. 



TABLE XVI: Terms in the EDF (30) that are of sixth order, 
depend on time-odd densities, and are built from the scalar 
nonlocal density p(r,r'). 



No. c; 



n L' V J 



Pn'L'v' J' 



D„ 



PnLvJ 



65 
66 
67 
68 
69 
70 
71 
72 



^1101 

'-^40,1101 



^1101 
(^20,3101 



Go' 



-,1101 

-'22,3101 
^1101 
'-'22,3303 



^00,5101 

^3101 

'-'00,3101 

^3303 

'^00,3303 



[Ml 


N^]o 


[Mil 


[Ml 


[vvio[vvi2 


[Mil 


[Mil 


[vvio 


[[fcfcloMli 


[Mil 


[vvi. 


[[fcfcloMli 


[Mil 


[vvi. 


[[k[kk],]^p]^ 


[Mil 


1 


[[fcfcloMli 


[[fcfeloMli 


1 


[[fcfcloMli 


[[k[kkU,p], 


1 


[[k[kk],],p]. 



TABLE XVII: Same as in Table XVI but for terms that are 
built from the vector nonlocal density s{r,r'). 



No. 


/~in' L' v' j' 


Pn'L'v'J' 


Dn,! 


PnLvJ 


73 


r'OOll 
'-'60,0011 


[s]. 


[vv];; 


[s]. 


74 


,-,0011 
'-'62,0011 


[^ll 


[VV]^[VVl2 


[^ll 


75 


^0011 

'-'40,2011 


[^ll 


[vvi^ 


[[fcfclo^li 


76 


^0011 

'-'42,2011 


[^ll 


[VV1JVV12 


[[fcfclo-sli 


77 


^0011 
'-'40,2211 


[*']l 


[vvif, 


[[kk],s]. 


78 


rMUl 
'-'42,2211 


[*']l 


[VV1JVV12 


[[fcfc]2^1i 


79 


/-,0011 

'-'42,2212 


[s]. 


[VV1JVV12 


[[kk],s]2 


80 


^0011 
'-'42,2213 


[s]j 


[VVIJVVI^ 


[[kk],s]^ 


81 


^0011 
'-'44,2213 


[^ll 


[V[V[VV]2]3]4 


[[fcfc]2«l3 


82 


^0011 

'-'20,4011 


[s]. 


[vvi„ 


[[kk]ls]. 


83 


^0011 

'-^22,4011 


[^ll 


[VV12 


[[kk]ls]. 


84 


^0011 
'-'20,4211 


[^ll 


[VVlo 


[[kk],[kk],s]^ 


85 


^0011 
'-'22,4211 


[s]. 


[VVl, 


l[kkUkk],sl 


86 


^0011 
'-'22,4212 


[^ll 


[VV12 


l[kk]olkk],s]2 


87 


^0011 
'-'22,4213 


[^ll 


[VV12 


[[kk],[kk],s]. 


88 


^0011 

'-'22,4413 


[s]. 


[VVl, 


mikm^i^is]^ 


89 


/-,0011 
'-'00,6011 


[*ll 


1 


imisi 


90 


^0011 
'^00,6211 


[s]j 


1 


[[kk]f,[kk],s], 



effective NN force from the bare NN force is the unitary 
correlation operator method (UCOM) [35]. The use of 
the UCOM scheme, however, leads to a nonlocal effective 
interaction even if the bare interaction would have been 
a local one. 

But rather than discussing to which extent gauge sym- 
metry is conserved or broken in nuclei we aim to provide 
a theoretical framework where different choices can be 
accommodated. Because several successful phenomeno- 
logical forces (e.g. Skyrme and Gogny [36]) are invariant 
under the gauge transformation, this symmetry consti- 
tutes a natural starting point in the search of improved 
EDFs. To which extent gauge symmetry is violated for 
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TABLE XVII: continued. 



No. 


^n' L'v' J' 
'^mI,nLvJ 


Pn'L'v'J' 


D^i 


pnLvJ 


91 


^2011 
^20,2011 


[fcfelo^li 


[VV]„ 


[[kk\A, 


92 


^2011 
(-^22,2011 


[kk\Ai 


[VV], 


[[kk],s\ 


93 


^2011 
"-^20,2211 


[kk\Ai 


[VV]o 


[[kk]A, 


94 


^2011 
(-^22,2211 


[kk]oS\i 


[VV]^ 


[[kk]A. 


95 


^2011 
(-^22,2212 


[kk\,s\i 


[VV], 


[\kk]A2 


96 


^2011 
■-^22,2213 


[kk]os\i 


[VV], 


[[kk]^. 


97 


^2011 
*-^00,4011 


[kk\as\i 


1 


[[kk]ls]. 


98 


^2011 
(-^00,4211 


[kk\A, 


1 


[[kkUkk]^^ 


99 


^2211 
^20,2211 


[kk]Ai 


[vv]o 


[[kk]A, 


100 


^2211 
(-^22,2211 


[kk]^s\ 


[VV]^ 


[[kk]^s]. 


101 


^2211 
L'22,2212 


[kk\Ai 


[VV], 


[[kk]A2 


102 


^2211 
^22,2213 


[kk]Ai 


[VV], 


[[kk],s]. 


103 


^2211 
*-^00,4011 


[kk]^s\ 


1 


[[kk]ls\ 


104 


^2211 
"-^00,4211 


[kk\Ai 


1 


[[kkUkk]^, 


105 


^2212 
^20,2212 


[kk]A2 


[VV]„ 


[[kk]A2 


106 


^2212 
(-^22,2212 


[kk\^s\^ 


[VV], 


[[kk\A^ 


107 


^2212 
L'22,2213 


[kk\A2 


[VV]^ 


[[kk\A^ 


108 


^2212 
^00,4212 


[kk]A2 


1 


[[kkUkk]^^ 


109 


^2213 
(-^20,2213 


[kk\A., 


[VV]„ 


[[kk]Ai 


110 


^2213 
"-^22,2213 


[kk\^s\. 


[VV]^ 


[[kk]Az 


111 


^2213 
*-^ 00, 4213 


[kk]^.^ 


1 


[[kk]o[kk]A^ 


112 


^2213 
*-^00,4413 


[kk\,sl. 


1 


[[k[k[kkU^\,sl, 



TABLE XVIII: Same as in Table XVI but for terms that 
are built from the scalar nonlocal density p{r,r') and vector 
nonlocal density s{r,r'). 



No. c: 



n' L V J 



Pn'L'v' J' 



D„ 



pnLvJ 



113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 



^1101 

^51,0011 

^1101 

(^31,2011 

^1101 

(--31,2211 

^1101 

^31,2212 

,-.1101 

(-^33,2212 

,-.1101 

^-33,2213 

,-.1101 

^11,4011 

,-.1101 

L'11,4211 

,-.1101 
(-'11,4212 



,-.3101 
(-'11,2011 

,-.3101 
(-'11,2211 

,-.3101 
(-'11,2212 

,-.3303 
(-'33 nnii 



/-.3303 
(-'11,2213 
,-.5101 
(-'11,0011 



[Ml 

[Ml 

[Mil 

[Mil 

[Mil 

[Mil 

[Mil 

[Mil 

[Mil 
[[fcfcloMli 
[[fcfcloMli 
[[fcfcloMli 
[[fcfcloMli 

[k[kk],],p\, 
[k[kk\,],pl, 
[k[kk].^]^p].:^ 
[[fcfcloMli 



[VV],^V 
[VV]„V 
[VVloV 
[VVIqV 
[V[VV]2]3 
[V[VV]2]3 

V 

V 

V 
[VV]„V 

V 

V 

V 
[V[VV],]3 

V 

V 

V 



[[kk]^^ 
[[kk]^^ 
[[kk]^^ 
[[kk]^^ 
[[kk]ls]^ 
[[kk]o[kk]A^ 
[[kk]o[kk]A^ 

[[kk]^^ 
[[kk],^^ 

[s]i 
[[kk]^^ 
[[kk].^^ 



TABLE XIX: Numbers of local densities pnLvj, Eq. (23), of 
different orders, which enter into the EDF up to N^LO. Num- 
bers of local densities constructed from the scalar p{r,r') or 
vector s{r,r') nonlocal densities, and numbers of time-even 
and time-odd local densities, are given separately. In Ta- 
bles III and IV these densities are marked with stars (*). 



order from p 


from s 


T-even 


T-odd total 





1 


1 


1 


1 


2 


1 


1 


3 


3 


1 


4 


2 


2 


4 


2 


4 


6 


3 


2 


6 


6 


2 


8 


4 


2 


5 


2 


5 


7 


5 


1 


4 


4 


1 


5 


6 


1 


2 


1 


2 


3 


total 


10 


25 


19 


16 


35 



effective renormalized interactions is a question that can 
be investigated by comparing models using preserved and 
broken symmetries (see e.g. Rcf. [37]). 



1. Local gauge transformations of the nonlocal densities 

The gauge-transformed nonlocal densities read [24, 25, 
36] 

p{r,r') = e^['^('')-'^(''')]p(r,r'), (34) 

s'{r,r') = e*['^M-'^(^')]s(r,r'). (35) 

Since the local gauge transformations form a U(l) group, 
invariance with respect to transformations that arc of the 
first-order in gauge angles, [1 -|- iG] p{r, r'), where 



G{r,r') 



(36) 



is enough to ensure full gauge invariance. By Taylor ex- 
panding the exponential functions in eq. 34 and 35 after 
they are inserted in the functional one may, of course, 
also prove this fact explicitly. 

One specific type of gauge transformation is the 
Galilean transformation, for which the gauge angles de- 
pend linearly on positions, i.e., G{r,r') = p ■ {r — r')/ti, 
and which corresponds to a transformation to a reference 
frame that moves with velocity p/m. For this transfor- 
mation, only first-order derivatives of G survive, which 
makes Galilean invariance less restrictive than the full 
gauge invariance. 



2. Local gauge transformations of the local densities 

Let indices /3, 7, . . . = 1, . . . , 35 label primary local den- 
sities PnLvJ (23) listed with stars (•) in Tables III and 
IV, which enter the EDF at N^LO, as shown in Eqs. (30) 



11 



and (31). Using this notation, the linearized gauge trans- 
formation of one of the local densities can be written as 



= pp{r) + {[if„LzG(r, r'K(r, r')] ^U,, 
= pp{r) + p''p{r), (37) 

where the first term is the untransformed local density 
and the second term is the part affected by the gauge 
transformation. 

As an illustration, let us begin by considering the sim- 
pler case of Galilean transformation, and look at the term 
with n — 2, where only two relative momentum operators 
k appear. Operator k can be written as kp-\-kQ, with the 
first term acting only on p^{r,r') and the second term 
acting only on G{r,r'). Then we have. 

K2L = [fcfcli = [fcpfcpli + 2 [fcpfccli + [fccfccli • (38) 

When this is inserted into the expression for p^ , the last 
term can be dropped since only the first-order derivatives 
of G{r,r') survive for the Galilean transformation, and 
the first term disappears when one takes the limit of r = 
r' since G{r,r) ~ 0. Thus in this case we obtain 

p^(r) ^ »{2[[fc,fcG],G(r,r')p.(r,r')],,,}^^^, 

= ^^c^'[{[fcpP.]J,},^,,(V0)(r)]^, (39) 
J' 

where the second equation results from the vector re- 
coupling (note that G{r,r') is a scalar and (V(/))(r) is a 
vector) and cj' are the ensuing numerical coefficients. 

This example illustrates the main features of the 
derivation, namely, (i) for the Galilean transformations 
only terms with first-order derivatives of G{r, r') occur 
in the final expression for p^, (ii) local densities appear- 
ing in the sum are of one order less in derivatives than 
the density being transformed, and (iii) the tensor or- 
der is preserved so that a local density is transformed 
into a sum of densities which can couple with a vector to 



the same tensor order. This leads to the ansatz for the 
Galilean transformation, 

p^ (r ) = * ^ c(/3, 7) K ir) (V0) (r)] ,, , (40) 



where c(/3, 7) are numerical coefficients. Similarly, for the 
full gauge transformation the corresponding ansatz reads 

P$i^) = * E Cm/(/3,7) K(^) [Drni<f\i {r)]j . (41) 

In both cases, the numerical coefficients can be found 
by using the method outlined above, combined with a 
repeated use of the 6j-symbols. 

However, instead of using this method, it turned out 
to be more efficient to proceed in another way. First, 
by using symbolic programming [38] , we constructed the 
transformed densities p^{r) explicitly in terms of deriva- 
tives of the density matrices and the gauge angle. Then, 
from the resulting expression the ansatz (40) or (41) was 
subtracted, which gave equations for the numerical coeffi- 
cients by requesting that these differences must be identi- 
cally equal to zero. Because these equations must hold for 
all density matrices and gauge angles, we could randomly 
assume arbitrary values for these quantities and their 
derivatives. In this way, all linearly independent equa- 
tions for the coefficients could be obtained and solved an- 
alytically, again by using symbolic programming. The so- 
lutions were then double-checked by using the full forms 
of the densities. 



3. Galilean or gauge invariant EDF 

A Galilean or gauge-invariant EDF is the one which 
does not change upon inserting Galilean or gauge- 
transformed densities (37) into the energy density of 
Eq. (31). Since terms quadratic in G{r,r') can be 
dropped, the condition for the Galilean or gauge invari- 
ance reads 



/ 



j3„V- 



'*" 2^ C'^/,7 



([pf (r) [D„,ip^{r)]j]^+ [pp{r) [i?„,p^(r)] ^^, 



= 0, 



(42) 



r 



where C^j ^ is a short-hand notation for the coupling 

constants G^f'^^j^'lj, and the sum runs over all the terms 
in the energy density. 

The task now is to group together all proportional 
terms in Eq. (42). In doing so, we do not aim at obtain- 
ing an invariant energy density but an invariant EDF and 
total energy. Therefore, after densities (40) or (41) are 



inserted into Eq. (42), all terms must be integrated by 
parts to obtain some standard form, where terms equal 
through integration by parts become identical. 

Finally, Eq. (42) can be transformed into a sum of 
independent terms using recoupling. In this expression 
each term is multiplied by a specific linear combination 
of coupling constants C'^^^if^j. The Galilean or gauge 
invariance of the EDF then means that these linear com- 
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binations must all vanish. This gives a set of linear equa- 
tions that must be fulfilled for an invariant EDF. On the 
one hand, if a given coupling constant appears in none 
of these linear equations, the corresponding term of the 
EDF is invariant on its own, and the corresponding cou- 
pling constant is not restricted by the Galilean or gauge 
symmetry. On the other hand, for some coupling con- 
stants the only solution can be the value of zero, and 
then the corresponding term cannot appear in the in- 
variant energy density. 

Among all the remaining coupling constants, we may 
always select a subset of those that we will call the de- 
pendent ones, and express them as linear combinations 
of the other ones, which we will call the independent 
ones. This procedure is highly non-unique and can be 
realized in very many different ways. However, when the 
dependent coupling constants in function of independent 
ones are inserted back into the energy density (31), lin- 
ear combinations of terms appearing at each independent 
coupling constant will all be invariant with respect to the 
Galilean or gauge transformations. 

Then, the energy density of Eq. (31) takes the form 



n{T) = 



A^ 



/-in'L'v'j' /-in'L'v'j' / \ 
^ml .nLvJ^ml .nLvJy )i 



(43) 



where the sum runs over indices that correspond to un- 
restricted and independent coupling constants, which we 
jointly call free coupling constants. For a term in Eq. (43) 
that corresponds to an unrestricted coupling constant, 
we have G';^f,'„li'j(r) = Tl^t^nH'Ar). i-c, one term in 
the energy density of Eq. (31) is Galilean or gauge in- 
variant. For a term in Eq. (43) that corresponds to an 
independent couphng constant, G^j ^jf^j{r) is equal to 
a specific linear combination of terms T from the original 
energy density (31). These linear combinations are listed 
in Appendix C. 

We performed the analysis along these lines for energy 
densities of orders 0, 2 , 4, and 6, and the obtained results 
are listed in Appendix C. Derivations were performed by 
using symbolic programming [38] and employed the tech- 
nique of forming linear equations by randomly assign- 
ing values to local densities and their derivatives, which 
we also used above. Numbers of linearly independent 
Galilean or gauge invariant terms are listed in Tables V 
and VI, and plotted in Fig. 1. 

It turns out that only at orders and 2, i.e., for the 
standard Skyrme functional, all Galilean invariant com- 
binations of terms are also gauge invariant. At orders 4 
and 6, there are only 6 gauge invariant terms available, 
while the numbers of Galilean invariant terms equal to 
15 and 26, respectively. This is much less than the total 
numbers of terms at these orders, which are equal to 45 
and 129, respectively. Altogether, at N'^LO we obtain the 
EDF parametrized in general by 188 coupling constants, 
and by 50 or 21 coupling constants if Galilean or gauge 
invariance is assumed. If isoscalar and isovector channels 
are included, all these numbers must be multiplied by a 



factor of two. 



IV. ENERGY DENSITY AT N^LO WITH 

CONSERVED SPHERICAL, SPACE-INVERSION, 

AND TIME-REVERSAL SYMMETRIES 

In this section, we apply the results obtained above to 
the simplest case of spherical even-even nuclei [28] , where 
one can assume that the spherical symmetry, along with 
the space inversion and time reversal, are simultaneously 
conserved symmetries. In this case, all primary densities 
PnLvj (23), which we listed in Tables III and IV, must 
have the form [39] : 

PnLvj{r)^ Rjj{r)pnLvj{\r\), (44) 



where 



i?Jj(r) = [r[r...,[rr]2,...,]j_2]j (45) 



is the Jth-ordcr, rank-J stretched coupled tensor built 
from the position vector r in exactly the same way as 
the derivative operators DnL of Table I are built from 
the derivative V in the spherical representation (9), and 
PnLvj{\i'\) is a scalar function depending only on the 
length |r| of the position vector r. 

Indeed, due to the generalized Cayley-Hamilton 
(GCH) theorem, a rank-J tensor function of a rank-fc 
tensor must be a linear combination of all independent 
rank-J tensors built from that rank-fc tensor, with scalar 
coefficients. In the GCH theorem, tensors that differ by 
scalar factors are not independent. In our case, only 
one independent rank-J function Rjj{r) can be built 
from the rank-1 tensor (position vector r), which gives 
Eq. (44). The spherical symmetry assumed here is essen- 
tial for this argument to work, because many more inde- 
pendent rank-J tensors can be built when other "mate- 
rial" tensors (like, e.g., the quadrupole deformation ten- 
sor) are available. 

The spherical form of PnLvj{f) (44) requires that the 
following selection rule is obeyed: 



P={-1) 



(46) 



where P = (—1)" is the space-inversion parity defined in 
Eq. (26). For the time-even densities {T = 1), selection 
rule (46) does not impose any new restriction on local 
densities built from p (r, r') {v = 0), see Table HI. On the 
other hand, for local densities built from s (r, r') {v — 1), 
see Table IV, only the densities with L — J are allowed. 
In Tables HI and IV, all densities allowed by the con- 
served spherical, space-inversion, and time-reversal sym- 
metries are marked with bullets (•). One can see that 
they correspond to quantum numbers LvJ being equal 
to 000 or 202 [for densities built from p{r^r')] and 111 
or 313 [for densities built from s {r^r')]. Then, it is easy 
to select all allowed terms in the energy density — in Ta- 
bles VII-XVIII and XXII these are also marked with bul- 
lets (•). Numbers of such terms are listed in Table XX 
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TABLE XX: Numbers of terms defined in Eq. (30) of different 
orders in the EDF up to N^LO, evaluated for the conserved 
spherical, space-inversion, and time-reversal symmetries. The 
last two columns give numbers of terms when the Galilean or 
gauge invariance is assumed, respectively, see Sec. IIIB. To 
take into account both isospin channels, the numbers of terms 
should be multiplied by a factor of two. 



order 


Total Galilean Gauge 





1 


1 1 


2 


4 


4 4 


4 


13 


9 3 


6 


32 


16 3 


N^LO 


50 


30 11 



together with those obtained by imposing, in addition, 
the Galilean or gauge invariance. 

All results for the EDF restricted by the spherical, 
space-inversion, and time-reversal symmetries can now 
be extracted from the general results presented in Sees. II 
and III and Appendices B and C. However, in the re- 
mainder of this section we give an example of how these 
results can be translated into those based on the Carte- 
sian representations of derivative operators (18)-(22). In- 
deed, in this representation, all non-zero densities can be 
defined as: 



those marked in Table IV. However, one should note 
that rank-2 densities R2ab and R^ab are not proportional 
to P2202 and ^4202, respectively, and the rank-3 density 
Jsabc is not proportional to P3313. This is so, because 
they are defined in terms of the derivative operators 
(18)-(22), where appropriate traces have not been sub- 
tracted out. Nevertheless, linear relations between den- 
sities (47)-(56) and their spherical-representation coun- 
terparts pnLvj can easily be worked out and will not be 
presented here. 

Note also that the scalar densities R2 and R4 can be 
expressed as the corresponding sums of the rank-2 den- 
sities i?2a6 and R^ab, and the vector density J^a as that 
of J^abc- However, based on the results obtained in the 
spherical representation, we know that they have to be 
treated separately to give separate terms in the energy 
density. 

Again, based on the results obtained in the spherical 
representation, we can write the N'^LO energy density 
as a sum of contributions from zero, second, fourth, and 
sixth orders: 



7i = TYo + ^2 + '^4 + ^67 



(58) 



where 



and 







Ro = P, 


(47) 






R2 = k^P, 


(48) 






R2ab = kakhp, 


(49) 






Ri = k^p, 


(50) 






Riab = k^kakhp, 


(51) 






Re = k^'p, 


(52) 


J la 


— 


[k X S)a, 


(53) 


Jsa 


= 


k^{k X S)a, 


(54) 


Jzabc 


~ 


kakb{k X s)c + kbkc{k x s)a 






+ 


kckaik X S)fc, 


(55) 


J5a 


= 


k'^ik X s)a, 


(56) 



Ho — CqqRqRo, (59) 



H2 — C2o-RoAi?o + CQ2R0R2 

+ C°ii?oV-Ji,+Co\J^ (60) 



Energy densities TYq and 7^2 correspond, of course, to 
the standard Skyrmc functional [24, 36] with Cqq — C*, 
C^o - C^' + \C\ C^o2 = C\ C?! = C^'\ and C^i = 
C'^^. At fourth and sixth orders, these energy densities 
read 



where 



fe — / ^ KaKai 



(57) 



and the Cartesian indices are defined as a,b,c = x, y, z. 
To lighten the notation, in these definitions we have omit- 
ted the arguments of local densities, r, and limits of 
r' — r. 

The six local densities (47)-(52) are the Cartesian ana- 
logues of densities marked in Table HI with bullets (•), 
and the four local densities (53)-(56) are analogues of 



Tii 



C^qRqA Ro 
Cq^RoR^ 



C22R0AR2 



22 
CQ2R2R: 



D22R0 2_^ '^a^bR2ab + Do2 /_^ R2abR2ab 



ab 



ab 



+ D\^Ji • V (V • Ji) 

+ C°ii?oA (V • Ji) + C?3i?o (V • Ja) 

+ Cl^R2{V-Ji)+Dl^Y.^^-bVaJib 



(61) 



ah 
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Hq — C^qRqIS. Rf) + C42i?oA i?2 



C24Ro^R4 + CqqRqRq 
C22R2^R2 + CQ4R2R4 
-D42-R0A 2_^ ^a^bR2ab + -D24-^0 2^ Va Vhi?4a6 

-^22-^2 2^ ^a^bR2ab + £^22 /_^ R2ab^R2ab 
ab ab 

-^22 2^ R2ab^a^cR2cb + -^04 /^ R2abRiab 
abc ab 

Cli Jl • A2 Ji + CI3 Ji . AJ3 
Co5 Jl • J5 + C'oS'^S • «^3 

Dl^Ji • AV (V • Jl) + i^ig Jl • V (V • J3) 

-^23 2^ Jla^b^ cJzabc + -Do3 2^ J ZabcJ Zabc 



abc 
^0 D_ A2 , 



a6c 



C^ii?oA^ (V • Jl) + C^gi^oA (V • J3) 
C°5i?o (V • J5) + C|ii?2A (V • Jl) 

C23i?2(V-J3) + Ci\i?4(V-Ji) 

D33R0 2_^ "^ a"^ b^ cJ S,abc + -D13 2^ -R2afc Vc J3afc 



a 6c 



'ibc 



-D3I 2^ R2ab^^aJlb + -^13 2^ R2ab^ aJ[ 



3b 



Dii 2_^ Riab'^aJlb 
ab 

Ei,J2R2abVaVbiV-Ji) 

ab 



(62) 



The energy densities above are given in terms of the 
coupling constants C^^, D^^^, E";^^, and F^'„. The 
indices correspond to orders of derivatives indicated in 
the same way as for the spherical-representation coupling 
constants C'^/n'^^j- Linear relations between both sets 
of coupling constants can easily be derived and arc not 
reported here. 



tional with density-independent coupling constants con- 
tains 188 terms, while functionals restricted by Galilean 
and gauge symmetries contain 50 and 21 terms, respec- 
tively. If both isoscalar and isovector channels are in- 
cluded, all these numbers must be multiplied by a factor 
of two 

At the present stage of searching for precise, 
spectroscopic-quality nuclear functionals, extensions be- 
yond the standard Skyrme NLO form arc mandatory, see 
the analysis in Ref . [40] . These may include richer density 
dependencies [41, 42], higher-order derivative terms, as 
constructed in the present study, terms of higher powers 
in densities, richer forms of functional dependence be- 
yond simple power expansions, and possibly many other 
modifications. 

Further studies of higher-order energy density func- 
tionals requires constructing appropriate codes to solve 
self-consistent equations. Although this is a complicated 
problem, various techniques have already been developed 
that can be used here. First of all, expressions for mean 
fields must be derived by using the standard methods 
presented, e.g., in Refs. [24, 25, 43]. Obviously, such 
mean fields will involve derivative operators up to sixth 
order, so the connection with the one-body Schrodinger 
equation, discussed, e.g., in Ref. [28], will be lost. Nev- 
ertheless, all basis-expansions methods can still be used 
and their implementation will not be basically different 
than it was done up to now at NLO. Work along these 
lines is now in progress. 

This work was supported in part by the Academy of 
Finland and University of Jyvaskyla within the FIDIPRO 
programme, by the Polish Ministry of Science under 
Contract No. N N202 328234, and by the UNEDF Sci- 
DAC Collaboration under the U.S. Department of En- 
ergy grant No. DE-FC02-07ER41457. 



APPENDIX A: SYMMETRIES 



V. CONCLUSIONS 

In the present study, we constructed nuclear energy 
density functionals in terms of derivatives of densities 
up to sixth order. This constitutes the next-to-next- 
to-next-to- leading order (N'^LO) expansion of the func- 
tional, whereby, in this scheme, the contact and standard 
Skyrme forces provide the zero-order (LO) and second- 
order (NLO) expansions, respectively. The higher-order 
terms were built to provide tools for testing convergence 
properties of methods based on energy density function- 
als, within the spirit of effective field theories. 

At N^LO, depending on several options of using the 
energy density functionals, the numbers of free coupling 
constants are as follows. If one would like to include 
the density dependence of all the coupling constants (the 
option, which is not advocated here), one would have 
to use 658 different terms in the functional. Full func- 



Within the HF approximation, symmetry properties 
of the interaction carry over to symmetry properties of 
the total HF energy. It means that whenever the inter- 
action is invariant with respect to a symmetry group, 
the total HF energy is invariant with respect to trans- 
forming densities (or density matrices, in general) by the 
same symmetry group. However, the well-known phe- 
nomenon of the spontaneous symmetry breaking [32] may 
render densities and energy density themselves not invari- 
ant with respect to the symmetry group in question. We 
must, therefore, consider energy densities for symmetry- 
breaking densities. 

This is best illustrated by the EDF derived for the 
Skyrme interaction [25, 28]. The standard derivation 
treats the time-reversal and isospin symmetries in a dif- 
ferent way than space symmetries (space inversion and 
other point symmetries or space rotation) [24, 25, 36]. 
Indeed, for the time reversal, the nonlocal densities are 
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first split into the timc-cvcn and time-odd parts as 

p{r,r') ^ p+{r,r')+p-{r,r'), (Al) 

s{r,r') = s+{r,r') + s-{r,r'), (A2) 

where 

P±ir,r') - i[p(r,r')±p^(r,r')], (A3) 

s±{r,r') = i[s(r,r')±5^(r,r')], (A4) 



appear in p{r) and s{r) in the energy density of Eq. (A8), 



such that 



p±ir,r') = ±p±{r,r'), 
s'^{r,r') = ±s±(r,r'). 



(A5) 
(A6) 



The superscript T here means that the nonlocal densities 
are calculated for the time-reversed many-body states. 

Then, in the derivation of the HF energy density, only 
squares of the time-even and time-odd densities appear, 
because, due to the time-reversal symmetry of the in- 
teraction, the cross terms do not contribute. As a con- 
sequence, the energy density itself is time-even. Simi- 
larly, for the isospin symmetry, the densities are first split 
into the isoscalar and isovector parts, for which no cross 
terms contribute, and the obtained energy density is an 
isoscalar. In what follows, we call this kind of deriva- 
tion 'derivation after separation of symmetries', which 
implies that the symmetry-breaking terms are absent in 
the energy density. 

The derivation after separation of symmetries can be il- 
lustrated by considering the simplest term of the Skyrme 
interaction - just the contact force. 



Vs{ri,r2) ^toS{ri 



r2 



(A7) 



for which the energy density reads (we neglected the 
isospin degree of freedom, so only one type of particles is 
considered), 

nsir) = itop'(r) - ^tos^ir). (A8) 

This energy density is invariant with respect to the time 
reversal of local densities, p^{r) = p{r) and s^ {r) = 
-sir). 

Here, the coupling constant multiplying the time-even 
density, C = 2^0; is not independent of the coupling 
constant multiplying the time-odd density, C^ ~ —\^o- 
This fact is not related to the time-reversal symmetry 
but results from the vanishing range of the contact force. 
Proper treatment of the finite range corrections render 
these two coupling constants independent of one another. 
[15, 16]. Irrespective of zero or finite range, the isovector 
and isoscalar coupling constants are also independent of 
one another [24, 36]. 

For space symmetries, the standard derivation pro- 
ceeds in another way, namely, the energy density is deter- 
mined directly for the broken-symmetry HF state. For 
the space-inversion symmetry, for example, this means 
that both parity-even and parity-odd densities, 

Pp=^,{r) - \ [p{r) ± pP{r)] = \ [p(r) ± p[-r)\ ikS>) 

sp^^.ir) = i [s{r) ± s^{r)] = i [sir) ± s{-r)lAlO) 



Pir) = Pp=+i{r)+Pp^_-i_{r), 



s{r) 



=+i(^) + «P=_i(^)- 



'P=+ 



(All) 

(A12) 



The superscript P here means that the nonlocal and lo- 
cal densities are calculated for the space-inversed many- 
body states. We call this kind of derivation 'derivation 
before separation of symmetries', which implies that the 
symmetry-breaking terms are then explicitly present in 
the energy density. If the symmetry is broken, which in 
the case of space inversion corresponds to Pp^_i{r) ^ 
or Sp^_^(r) 7^ 0, then the energy density is noi invariant 
with respect to space inversion. 

The total energy, i.e., the integral of the energy density 
(1) is, of course, invariant with respect to space inversion, 
because the integration then picks up only the space- 
inversion-invariant parts of the integrand. Therefore, the 
energy densities derived before and after separation of 
symmetries are not equal, but they are equivalent. In the 
case of the space-inversion symmetry, the energy density 
(A8) derived after separation of symmetries reads 









(A13) 



This energy density is invariant with respect to space in- 
version and the coupling constants multiplying densities 
of opposite parities are not independent of one another 
(see also discussion in Ref. [30]). 

The same principle applies to other broken spatial sym- 
metries. For example, for the broken rotational sym- 
metry, the density can be split into the sum of terms 
belonging to different irreducible representations of the 
rotational group, which in this case corresponds to the 
standard multipole series [44], 



Pir) =^PAA<(r), 



for 



Px^{r)=p^^{r)Y^^{0,<P), 
where the multipole densities pXfi{r), 

px^{r)^ f dddMr)Yx^i9,4 



(A14) 



(A15) 



(A16) 



depend only on the radial coordinate r. Then, the first 
term in the energy density (AS), which is derived before 
separation of rotational symmetries, and which is not ro- 
tationally invariant, is equivalent to the following energy 
density derived after separation of rotational symmetries, 

^p(^) = ¥o Y. V2XTT[px{r)px{r)]o, (A17) 



(see Ref. [45] for an example application of this series). 
This energy density is rotationally invariant and the cou- 
pling constants multiplying different multipole densities 
are again not independent of one another. 
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We have presented a detailed analysis of the problem 
to arm ourselves with proper tools for discussing con- 
struction of EDFs in situations where there is no un- 
derlying interaction. Then, the only consideration is the 
requirement of invariance of the total energy with respect 
to all symmetries usually conserved by nuclear interac- 
tions. We proceed with such a construction in two dif- 
ferent ways as described below. 



primary densities PnLvjix) (23) are equal to (—1)"+", see 
Eq. (25) and columns denoted by T in Tables III and IV. 
Similarly, T-parities of secondary densities Pmi.nLvj,Q{i') 
(24) are also equal to (—1)"+". Construction of the T- 
invariant energy density (A18) can now be realized by 
multiplying densities that have identical T-parities. 



2. Symmetry-covariant energy density 



1. Symmetry-invariant energy density 

Based on the derivation after separation of symme- 
tries, which we introduced above, it is clear that we can 
proceed by separating densities into irreducible represen- 
tations of all required symmetries and then building the 
energy density by taking scalar products separately in 
each of the representations. Such a construction gives a 
symmetry-invariant energy density. 



H^(r)=H(r), 



(A18) 



where Ti.^{r) denotes the energy density calculated for a 
many-body state transformed by the symmetry operator 
S. This guarantees the invariance of the EDF and to- 
tal energy (1) with respect to all considered symmetries. 
Such a strategy would also allow for using arbitrary, un- 
related to one another coupling constants in each of the 
irreducible representations. 

However, in practical applications, such a strategy was 
never up to now fully implemented — neither at N'^LO, for 
which the present study is the first attempt in the liter- 
ature, nor at NLO, which corresponds to the standard 
Skyrme functionals (see Sec. Ill A below). Only the time 
reversal and isospin symmetries were up to now treated 
in this way, and below we are going to follow the same 
path. 

For the time reversal, all local densities discussed in 
Sec. II C are either time-even or time-odd. Indeed, this 
simply follows from the facts [25] that 



p^{r,r') = p*{r,r') = 



\ry) 



s*(r, r') 



p{r',r), 
s{r',r), 



(A19) 



which give the time-even and time-odd parts in Eqs. (A3) 
and (A4) as 



P+{r,r') = 

p-{r,r') = 

s+{r,r') = 

s^{r,r') = 



p*+{ry) = 

p*_{r,r') = 

s*^{r,r') = 

s*_{r,r') = 



P+ir',r), 
P-ir',r), 
s+{r',r), 
s^{r' , r), 



(A20) 



i.e., pj^{r,r') and s^{r^r') are real symmetric functions 
and p- (r, r') and s+(r, r') are imaginary antisymmetric 
functions. Moreover, the relative momentum operator k 
(6), which defines derivative operators KnL, is imaginary 
and antisymmetric with respect to exchanging variables 
r and r' . Altogether, it is easy to see that T-parities of 



Treatment of space symmetries in the construction of 
EDFs is another kettle of fish. Here, we base our consid- 
erations on the derivation before separation of symme- 
tries, which we introduced above, and on the fact that 
invariance of the energy density itself is not a prerequi- 
site for the invariance of the EDF. In fact, the EDF and 
total energy (1) are invariant with respect to symmetry 
S also when the energy density is covariant with S, i.e., 



n^{r) = n{S+rS), 



(A21) 



where S'^rS denotes the space point transformed by sym- 
metry S (see also discussion in Ref. [32]). Indeed, due to 
the fact that space integrals are invariant, we have 

f d^rn{S+rS) ^ f d^rn{r), (A22) 

which guarantees invariance of the EDF and total energy. 
For the space-inversion symmetry, we have 



n^{r) = W[p^(ra,rV)] 
= W [p(-rcr, -rV)] 
n{P+rP) = H{-r), 

and the covariance condition (A21) reads 

W [p(-rcr, -r'cr')] =n{-r). 



(A23) 
(A24) 



(A25) 



It is now essential to realize that the arguments of the 
density matrix, p{—ra, —r'a') on which the energy den- 
sities in Eq. (A25) depends, are the same on both sides 
of Eq. (A25). The covariance condition then tests only 
the parity of all other operators that may appear in the 
definition of local densities. Therefore, to each primary 
density PnLvjif) (23) we may attribute P-parity corre- 
sponding to the P-parity of the operator KnL only, which 
is equal to (—1)", see Eq. (26) and columns denoted by 
P in Tables III and IV. This attribution is performed re- 
gardless of space-inversion properties of the nonlocal den- 
sities, i.e., regardless of whether the parity of the many- 
body state is conserved or broken. Similarly, P-parities 
of secondary densities Pmi,nLvj,Qif') (24) are equal to 
(—1)"+™. Construction of the P-covariant energy den- 
sity (A21) can now be realized by multiplying densities 
that have identical P-parities. 

Construction of a rotationally covariant energy density 
can be performed in an entirely analogous way. We must 
only ensure, that all tensor operators used in constructing 



all terms of the energy density are always coupled to total 
angular momentum (rank) zero. This coupling proceeds 
regardless of any transformation properties of nonlocal 
densities with respect to rotation, because again, their 
rotated space arguments appear on both sides of the co- 
variance condition (A21). 

It is obvious that this is the correct procedure to follow 
when the rotational symmetry is not broken, and nonlo- 
cal densities p{r,r') and s{r,r') are scalar and vector 
functions of their arguments, respectively. In fact, this 
is how we refer to these densities throughout the entire 
paper, seemingly forgetting that the rotational symme- 
try can be broken, and that these functions can then 
have no good tensor properties with respect to rotation. 
Nevertheless, in view of the covariance condition (A21), 



J 
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these rotational properties of broken-symmetry nonlocal 
densities are irrelevant for the construction of the energy 
density. 



APPENDIX B: PHASE CONVENTIONS 

In the present study, we use four elementary building 
blocks to construct the EDF, namely, the scalar and vec- 
tor nonlocal densities, p{r,r') and s{r,r'), along with 
the total derivative V and relative momentum k (6). 
Spherical representations of the building blocks can be 
defined by using standard convention of spherical tensors 
[26] as 



Poo{r,r') 
'^l,^l={-lfi,l} 



PpP{r,r'), (Bl) 

Ps { ^ (s,(r, r') - isyir, r')) , s,(r, r'), ^ (s,(r, r') + is,(r, r'))} , (B2) 

P^{j^i^^-^^y).^z,^{^.+i^y)}, (B3) 

Pk{i^i'^--^'^y)^'^-^^i'^^ + ^'^v)}> (B4) 
I 



where pp, ps, p-^, and p^, are arbitrary phase factors, 
\pp\ — \ps\ — \py\ = \Pf.\ = 1. These phase factors define 
the phase convention of the building blocks, and can be 
used to achieve specific phase properties of densities and 
terms in the EDF, as discussed in this Appendix. 

In order to motivate the best suitable choice of the 
phase convention, in Tables XXI and XXII we present 
relations between the spherical and Cartesian represen- 
tations of densities and terms in the EDF, respectively. 
All NLO densities in the Cartesian representation, which 
are listed in Table XXI, are real. It is then clear that the 
phase convention, which would render all NLO densities 
in the spherical representation real does not exist. How- 
ever, for the phase factors Pp, Ps, p-^, and p^, equal to ±1 
or ±i, in the spherical representation all NLO densities 
and terms in the EDF are either real or imaginary. 

Among many options of choosing the phase convention, 
in the present study we set 



Pp = +1, Ps = -i, Pv 



and Pi^ = 



(B5) 



This choice is unique in the fact that all scalar densi- 
ties and all terms in the EDF arc then characterized by 
phase factors +1 connecting the spherical and Cartesian 
representations, see the last columns in Tables XXI and 
XXII. This allows for the closest possible relationships 
between both representations, which may facilitate the 
use of the spherical representation as it is introduced in 
the present study. In particular, relations between cou- 
pling constants up to NLO (Table XXII) and standard 



coupling constants in the Cartesian representation [24] 
then read, for terms depending on time-even densities: 



^0000 
•-^OOjOOOO 

/-.OOOO 
'-^20,0000 

/-.OOOO 
^00,2000 

/-lino 

*-^00,1110 

/-.nil 

*-^00,llll 

/-.1112 
'-^00,1112 

/-.OOOO 

^11,1111 



^/3C^ 



jcn 



and for terms depending on time-odd densities: 



/-.OOll 
"-^00,0011 

^1101 
^00,1101 

^0011 
^20,0011 

^0011 
^22,0011 

^0011 
■-^00,2011 

^0011 
'-^00,2211 

^1101 



^/3C^ 

3 (c^^ + ic^) 4 



c^ - c 



Vs 



(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 
(B14) 
(B15) 
(B16) 
(BIT) 
(B18) 
(B19) 



At the same time, all vector densities in Table XXI and 
vector operators in Eqs. (B2)-(B4) are consistently char- 
acterized by phase factors —i connecting the spherical 
and Cartesian representations. 
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TABLE XXI: Spherical and Cartesian representations of local densities (24) up to NLO. Only scalar densities and the /i = 
components of vector densities are shown. Numbers in the first column refer to numbers of primary densities (23) shown in 
Tables III and IV. The last column shows factors preceding densities in the Cartesian representation evaluated for the phase 
conventions of Eq. (B5). Time-even densities are marked by using the bold- face font. 



No. 


Pm.I,nLvJ,Qfi 










Cartesian representation [24, 25, 30] 


Phase 


1 


Poo, 0000, 00 


= 


[p]oo 


= 


Pp 


P 




+ 1 




Pll,0000,10 


= 


[Vp]io 


= 


PvPp 


\7zp 




— i 




P20,0000,00 


= 


[[VV]op]oo 


= 


- P%Pp 


Ts^P 




+ 1 


2 


poo, 1101, 10 


= 


[kp]io 


= 


PkPp 


jz 




— i 




pll,1101,00 


= 


[V [ftp] 1 loo 


= 


- PvPkPp 


73^-i 




+ 1 




Pll,1101,10 


= 


[V[Mi]io 


= 


iPvPkPp 


T^(Vxi). 




—i 


3 


Poo, 2000, 00 


= 


[[fcfe]op]oo 


= 


-pIpp 


73 (--iM 




+ 1 


17 


poo, 0011, 10 


= 


[s]io 


= 


Ps 


Sz 




— i 




Pll,0011,00 


= 


[Vs]oo 


= 


- PvP= 


73^- 




+ 1 




Pll,0011,10 


= 


[Vs]io 


= 


iPvPs 


;^(^x-). 




—i 




P20,0011,10 


= 


[[VV]os]io 


= 


' P^Ps 


T^^«^ 




—i 




p22,0011,10 


= 


[[VVl2S]l0 


= 


- P%Ps 


^(3V.V.s- 


As.) 


—i 


18 


Poo, 1110, 00 


= 


[fcs]oo 


= 


- PkPs 


^,^''' 




+ 1 




Pll,1110,10 


= 


[V[fcs]o]io 


= 


- PvPkP^ 


;^v./°) 




—i 


19 


Poo, 1111, 10 


= 


[ks] 10 


= 


iPkPs 


^.J^ 




—i 




Pll,llll,00 


= 


[V[fcs]i]oo 


= 


-iPvPkPs 


ify--^ 




+ 1 




Pll,llll,10 


= 


[V[Mi]io 


= 


- PvPkPo 


\ (V X J)^ 




—i 


21 


poo, 2011, 10 


= 


[[A;fc]os]io 


= 


~pIp- 


^AT^-\^^^) 




~i 


22 


poo, 2211, 10 


= 


[[A;fe]2s]i() 


= 


-pIp^ 


-h, (3^^ - 1^^ 


V-s-T. + iAs.) 


—i 



TABLE XXII: Spherical and Cartesian representations of terms in the EDF (30) up to NLO. The last column shows factors 
preceding terms in the Cartesian representation evaluated using the phase conventions of Eq. (B5). Integration by parts 
was used to transform s ■ W • s into —(V • s)^, which is the term used previously in Refs. [24, 30]. Coupling constants 
corresponding to terms that depend on time-even densities are marked by using the bold-face font. Bullets (•) mark coupling 
constants corresponding to terms that do not vanish for conserved spherical, space-inversion, and time-reversal symmetries, see 
Sec. IV. 



No. 


r^n' L' v' j' 
^ml ,nLvJ 


[Pn'L'v'J'lDmlpnLvj] j, 


]o 


Cartesian representation 


Phase 


1 


. /-.OOOO 
• '-'00,0000 


Mo 


= Pp 


p' 


+1 


2 


^0011 
"-^ 00, 0011 


[ss]o 


= -ri 


V3^ 


+1 


3 


, ^0000 
• '-'20,0000 


[p[[VV]„p]„]„ 


= - PvPp 


-kP^P 


+1 


4 


, ,-,0000 
• '-'00,2000 


[p[[kk]oP]a]o 


2 2 

= - PfcPp 


^Ap--\p^p) 


+1 


5 


/-.lllO 
<-^00,1110 


[[ksUks\X 


2 2 
= PkPs 


\{J'''Y 


+1 


6 


^1111 

• '-'00,1111 


[[MJMJo 


2 2 
= PkPs 


^^^ 


+1 


7 


/-,1112 
•-'00,1112 


[[ksUks]^], 


2 2 
= PkPs 


1 Y^ r(2) ^(2) 

^ Z^ab ■^ab ■^ ab 


+1 


8 


_ ,-,0000 

• '-'11,1111 


[p[V[Mi]o]o 


= -i-PvPkPoPP 


ieP^-J 


+1 


9 


^1101 
(-'00,1101 


[MJMilo 


2 2 
= - PfePp 


71^-^ 


+1 


10 


^0011 
<-'20,0011 


[s[[VV]„s]Jo 


2 2 
= P^Ps 


IsAs 


+1 


11 


^0011 
'-'22,0011 


[s[[VV]2s]Jo 


= pIp's 


-i((V.s)^ + isAs) 


+1 


12 


^0011 
'-'00,2011 


[s[[fcfc](,s]J„ 


= Plpi 


i(s.T-isAs) 


+1 


13 


^0011 
'-'00,2211 


[s[[kk],s],]^ 


2 2 
= PkPs 


-Is-T+^^sAs) 


+1 


14 


^1101 
'-'11,0011 


[[Mi[Vs]j„ 


= -iPvPkP^Pp 


i?.3-^x^ 


+1 
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Phase conventions (B5) also lead to very simple phase 
properties, which our spherical tensors have with respect 
to complex conjugation. Indeed, spherical tensors (Bl)- 
(B4) obey standard transformation rules under complex 
conjugation [26], 



AX 



Pa(-1)^-^Aa 



(B20) 



where Pa — ±1- For nonlocal densities (B1)-(B2), 
Eq. (B20) holds separately for their time-even and time- 
odd parts, split as in Eqs. (Al) and (A2). Using 
Eqs. (A20) we then have 






-Pp 

-pI 

„2 



Pp- = -P% 

Ps- = -pI 



(B21) 



which for the phase convention of Eq. (B5) reads 



P. 



p+ 



-1 , P. 



Ps, = -1 



p- ^ 

p,_ = 



(B22) 



P^ 



-1 , Pk 



-1. 



Standard rule (B20) propagates through the angular 
momentum coupling, i.e., if signs Pa and Pa' character- 
ize tensors Ax and A'y, respectively, then the coupled 
tensor, 

A'i„^„ = [A,A',,]y,,„ = J2 <aV^a^^'vp', (B23) 

Pfi' 

is characterized by the product of signs Pa" = PaPa' ■ 
Therefore, coupled higher-order densities (24) are char- 
acterized by signs, 



P, 



PmI,nL'uJ,Q 



P^P^P.T. 



(B24) 



where w = or 1 denotes the scalar or vector density, p or 
s, respectively, and T = -|-lorT = — 1 denotes the time- 
even or time-odd density. However, symmetry conditions 
(A20) require that powers of the k derivative determine 
the time-reversal symmetry of each local density, so that 
T = (-1)"+". From Eqs. (B21) we then obtain 



^.™.,„.„.,<3 = (-l)'"+"+XVfe>^. (B25) 

which for the phase convention of Eq. (B5) reads 



^PmI,nL^.J,Q ^ '^ 



(B26) 



for all densities. Therefore, the phase convention of 
Eq. (B5) ensures that scalar densities and all terms in 
the EDF are always real. 



APPENDIX C: RESULTS FOR THE GALILEAN 

OR GAUGE INVARIANT ENERGY DENSITY 

FUNCTIONAL 

As discussed in Sec. IIIB 3, when the Galilean or gauge 
invariance is imposed on the EDF, this induces specific 



TABLE XXIIL Numbers of unrestricted, vanishing, indepen- 
dent, and dependent coupling constants in the EDF at zero, 
second, fourth, and sixth orders. Left and right columns cor- 
respond to the Galilean and gauge symmetries imposed, re- 
spectively. 







Galilean 






Gang 


e 




order 





2 


4 


6 





2 


4 


6 


unrestricted 


2 


3 


3 


3 


2 


3 


3 


3 


vanishing 




















27 


100 


independent 





4 


12 


23 





4 


3 


3 


dependent 





5 


30 


103 





5 


12 


23 



constraints on the coupling constants and terms of the 
functional. We pointed out that there can be three dis- 
connected classes of terms in the EDF with related prop- 
erties of the coupling constants: 

1. Terms that are invariant with respect to the 
Galilean or gauge transformation, and, therefore, 
the corresponding coupling constants are not re- 
stricted by the imposed symmetries. 

2. Terms that cannot appear in the energy density 
when the Galilean or gauge symmetry is imposed, 
and, therefore, the corresponding coupling con- 
stants must be equal to zero. 

3. Terms that can appear in the energy density only 
in certain specific linear combinations with other 
terms. This means that the coupling constants cor- 
responding to these terms must obey specific linear 
conditions. We then distinguish: 

(a) independent coupling constants, which mul- 
tiply invariant combinations of terms and, 
therefore, their values are not restricted by the 
imposed symmetries and 

(b) dependent coupling constants, which are equal 
to specific linear combinations of independent 
coupling constants, and, therefore, their val- 
ues are in this way restricted by the imposed 
symmetries. 

Division into the sets of independent and dependent 
coupling constants is not unique, and below, in each 
case, we present only one specific choice thereof. 

In Table XXIII we show numbers of unrestricted, vanish- 
ing, independent, and dependent coupling constants that 
appear at a given order when cither Galilean or gauge 
symmetry is imposed. In what follows, we use the name 
of a free coupling constant to denote either the unre- 
stricted or independent one. Indeed, in the Galilean or 
gauge invariant energy density (43), these two groups of 
coupling constants become free parameters. 
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Below we simultaneously discuss the Galilean and 
gauge symmetries. In doing so, we use the fact that the 
Galilean symmetry is a special case of the gauge sym- 
metry, and, therefore, the latter may impose more re- 
strictions on the EDF than the former. At NLO, this is 
not the case, and the Galilean and gauge symmetries im- 
pose, in fact, identical restrictions on the EDF [25, 36]. 
However, at higher orders, restrictions imposed by the 
Galilean and gauge symmetries are very different. 

Both zero-order terms in the EDF, which correspond to 
the contact interaction, are Galilean and gauge invariant, 
i.e., these symmetries do not restrict the form of the EDF 
at LO. In the three following sections we give results for 
second, fourth, and sixth orders, respectively. 



Again, terms that depend on time-even densities are 
marked by using the bold-face font. Altogether, 7 free 
coupling constants (3 unrestricted and 4 independent) 
define the Galilean or gauge invariant EDF at second or- 
der, cf. Table VI. 



Fourth order 



1. Second order 

At second order, we obtain the same restrictions of 
the EDF as those already identified for the Skyrme func- 
tional, see Ref. [24] for a complete list thereof. Then, 5 
dependent coupling constants are equal to specific linear 
combinations of 4 independent ones: 
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These relations are obtained by imposing either Galilean 
or gauge invariance. In Eqs. (C1)-(C5), coupling con- 
stants corresponding to terms that depend on time-even 
densities are marked by using the bold- face font. The 
same convention also applies below. 

At this order, the Galilean or gauge invariant energy 
density of Eq. (43) is composed of three terms corre- 
sponding to unrestricted coupling constants: 
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and of four terms corresponding to the independent cou- 
pling constants: 
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At fourth order, imposing either Galilean or gauge 
symmetry forces 12 dependent coupling constants to be 
specific linear combinations of 3 independent ones: 
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At this order, there are 3 stand-alone Galilean and gauge 
invariant terms: 
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and 3 Galilean and gauge invariant linear combinations 
of terms, corresponding to the 3 independent coupling 
constants: 
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Altogether, these 6 free coupHng constants (3 unre- 
stricted and 3 independent) occur in both the Gahlean 
and gauge invariant EDF at fourth order, cf. Table VI. 

Apart from these 6 free and 12 dependent coupling 
constants, the gauge invariance requires that all the re- 
maining 27 coupling constants are equal to zero. These 
27 constants are allowed to be non-zero if the Galilean 
symmetry is imposed instead of the full gauge invariance. 
Then, there are 18 dependent coupling constants that are 
forced to be linear combinations of 9 independent ones: 



c 



^0000 
'^•>n ■?.ni 



^1101 



c. 



0000 
22,2202 



^1110 
*-'20,1110 

^1110 
•-^22,1112 



--rllll 
^20,1111 



--<1111 
^22,1111 



1101 
22.1101' 



1/-.0011 



3^20,2011 3 V O<-'20,2211; 



^1111 2 --flll2 

•-^22,1112 yf^22,1112 



_0 /iS^'OOll 
Z,\/ — L'22,2213) 



5^0011 
3 '-'20,2211 



1 ^0011 
"yf "-'20,2011' 



a 



0011 
22,2011 



_ 1 /^llll 
y3«-22,1112 



5/^0011 
L/99 99- 



3 V 0^-^20,2011 

2/-rllll _ 2 

3*-'22,1112 77 

5 r'^011 

21^22,2213' 



3/^1112 
7^22,1112 



1^0011 
" 6^20,2211' 



(C31) 
(C32) 
(C33) 

(C34) 
(C35) 

(C36) 

(C37) 

(C38) 



0011 
22,2211 



X-.0011 
"-^22,2212 



/-<0000 
^%^ 11 



/-<00C 
*-^ll, 



3111 



/^2000 
'-^11,1111 

>^2202 
'-^11,1111 

^2202 
'-'11,1112 

^1101 
"-^11,2011 

X-.1101 
"-^11,2211 

/-(3101 
"-'11,0011 



1 /r^llll 
6 V 0»->22,1112 

I 8 r^OOll 
"1":^^ "-'22, 2213' 

1/^^1111 
2VOL>22,1112 

,0 /s'r'OOll 
"T^Y 7"-^22,2213 

/-fllOl 
"-'31,0011' 

_ /J^llOl 
Y 5"-'ll,2212' 



5/^1112 
7*-^22,1112 



3/-«1112 
7*-^22,1112 



'5^*1101 
l"-'ll 

LOl 
73 -'11,2212' 



1 ^11 

73^11 

-illOl 



3^11,2212' 
1101 



"-'11,2212' 



- Vs'-'ii 



,2212' 



1 ^1101 
75 "-11,2212' 



3r'iioi 
"-'11 



2212- 



(C39) 

(C40) 
(C41) 

(C42) 

(C43) 

(C44) 
(C45) 

(C46) 

(C47) 

(C48) 



Finally, we list 9 combinations of terms that are in- 
variant with respect to the Galilean symmetry and cor- 
respond to the independent coupling constants: 
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3. Sixth order 

An entirely analogous pattern of terms and cou- 
pling constants appears at sixth order. Imposing either 
Galilean or gauge symmetry forces 23 dependent cou- 
pling constants to be specific linear combinations of 3 
independent ones:. 
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At this order, there are 3 stand-alone Galilean and gauge 
invariant terms: 
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and 3 Galilean and gauge invariant linear combinations 
of terms, corresponding to the 3 independent coupling 
constants: 
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Altogether, 6 free coupling constants (3 unrestricted and 
3 independent) occur in both the Galilean and gauge in- 
variant EDF at sixth order, cf. Table VI. 

Apart from the 6 free and 23 dependent coupling con- 
stants, at sixth order the gauge invariance requires that 
all the remaining 100 coupling constants are equal to 
zero. These 100 constants are allowed to be non-zero 
if the Galilean symmetry is imposed instead of the full 
gauge invariance. Then, there are 80 dependent coupling 
constants that are forced to be linear combinations of 20 
independent ones: 
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L/9 



(_/9 



2213 



4 /^2011 
5^5 '-"22,2212 

^5V I4C22.3? 



c; 



X-.2213 
"-^22,2213 



/-<0000 
'-"51,1111 

/-.OGOG 
*-^ai .qi 1 1 



33,3313 
0000 



y-fOOOO 

'-"11,5111 

--.2000 
'-"31,1111 

/-.2000 
'-"11,3111 

/-.2202 
'-"31,1111 

,-.2202 
'-33,1111 

,-.2202 
'-31,1112 

,-.2202 
'-33,1112 



1 ,-.1111 
/15 '-22,3112 



22,33131 



(C130) 



=c. 



2011 



25V^ 22,2212 - 

I 4 r^llll I 

^ 5V2T '-'22,3112 + 

2 /-(2011 1 



3 ,-1110 
7 '-22,3112 



2 ,-.1111 
5 '-22,3313' 



;(C131) 



/15 



.c 



20.2011 



2^3 



C20.22111 (C132) 



cV (09929 



ir 



nil 

22,3112 



3 /-.nil 

TTo^' 



22,33131 

2 ,-.2011 
35 "-22,2212 



(C133) 



14,-.llll 
Is "-22,3112 



3/-I1111 
5'-22, 33131 



3 "-20,2011 



(C134) 



/105 



^20,2211' (C135) 



6 ,-.2011 
7'-22,2212 



1 ,-.1110 
7l4 '-22,3112 



I 2 /2,-.llll 
' T5 V 7'-22,3112 

,-.1101 
"-51.0011' 



3^1111 



3r'1101 

5 "-31,2212) 



C 



1101 
33,2213' 

J_ 
14 



1 c/^3101 
^^"-11.2212' 



S^illOl 

3 "-31,2212' 



_ /5/-<3il 



,2212' 



1 X-.1101 

yg "-31,2212' 

V "J"-33,2213' 

,-,1101 
"-31,2212' 

V •J"-33,2213' 



22,331(3^136) 

(C137) 
(C138) 
(C139) 
(C140) 

(C141) 

(C142) 
(C143) 
(C144) 
(C145) 
(C146) 
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'-^11,3111 

^2202 
^11,3112 



O- 



3101 
31,0011 



^3101 
'-'11,2011 



C{ 



3101 
1,2211 



^3303 
•-"33,0011 

,^3303 
'-'11,2212 

/-.3303 
'-"11,2213 

^5101 
•-"ILOOll 



1 ^3101 
^3 "-"11,2212' 

^3101 
'-"11,2212' 



•-^11,3312 — 



^11,3313 

>-,4000 

*-^ll,llll 

/^4202 
^11,1111 

/-f4202 
*-^ll,1112 

^1101 
"-"31,2011 

^1101 
'-'31,2211 

^1101 
'-"33,2212 

^1101 
'-"11,4011 

^1101 
'-"11,4211 

^1101 
'-"11.4212 



b x^aiOi 

21^11,2212; 

10^3101 
2T "-^11, 2212; 



^31C 
"-•ll. 



22121 



/-<31C 

f'-"ii.: 



7V3 ""11.2212; 

Sr'SlOl 

7"-"ll,2212; 

_ /5"/^1101 
Y 3'-"31,2212; 

1 ^1101 
:^ '-'31, 2212) 

_2./9riioi 

zv ^'-^33,2213' 



r<310 
'-'11.2 



5 /^3101 

7:73*-"ii,: 



3^11,2212! 

773 -11, 2212 > 

5^3101 

7'-^ll,2212) 



3/-11101 

5 '-"31,22121 



_ 7 5/^3101 
y 3'-"ll,2212) 

1 ,^3101 
7f '-'11,2212> 



V3 

^1101 

'-"33,2213' 



2 / 5 ^3101 

3 V 2T'-^11,2212' 

_2 
3' 



2 710^3101 

21 "-"11, 2212 I 



-I r^3101 

1^ "J '-"11, 2212- 



(C147) 
(C148) 

(C149) 
(C150) 

(C151) 
(C152) 
(C153) 
(C154) 
(C155) 
(C156) 
(C157) 
(C158) 
(C159) 

(C160) 

(C161) 
(C162) 
(C163) 
(C164) 

(C165) 
(C166) 



Finally, we list 20 combinations of terms that are in- 
variant with respect to the Galilean symmetry and cor- 
respond to the independent coupling constants: 
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G: 



2202 
20,2202 



G 



1110 
42.1112 



^1110 
^22,3112 



^1111 
^42,1111 



G 



1111 

22,3112 



K/cT.2000 
2^0-1 20,2000 



Tn2202 
-'^20,2202 



6 rrnllOl 

"^-'20,3101 



TnlllO I 1 rplll2 _ IthOOII 

-'■42,1112 + :^-'42, 1112 2 42,2011 

rplll2 I TnlllO I 2 ^1112 

-^22,3110 ' -'22,3112 ' :^-'22,3112 " 



3 T-iOOOO 
^^-^ 20,40001 

1 ^0011 3 

"^-^42,2211 2 



-illlG 
21-^22,3312 



3 7^0011 
35 -* 42,2213' 



-t1112 
3-^22,3312 



(C167) 
(C168) 



1 V 0-t22,3 



2 
3314 



5 7^2011 
"T2 22,2011 



1,/Ft2011 _ IrrnOOll _ l^/^T^OOll 3 /ISt^OOH 

gV 0-^22,2211 4-'^22,4011 7VO-t22,4211 TIVT 22,4213 



1 ,/Fti0011 _ 1^^2211 
^V 0^22,4413 3 -'22,2211 



2 

pllll 

' 42, 1111 



'15ti2011 
T"-' 22,2213 



3ti2211 
7-^22,2213 



zv •3-142,1112 



1 7^2213 
"^-'22,2213' 

nOOll 



^-.1112 I 4 rrnOOll 4 rpOOll _6 710011 

7-^42, 1112 + :^-'42, 2011 :^%-'42,2211 "^^ff-* 42,2213 J 



1 T-.1111 
^-t 22,3111 



1 rril 



'^22,3111^-^22,3112 "^-^22,3112 



-illlO 

51-^22,3312 



3-^22,3312 



(C169) 
(C170) 



^1111 

'-^22,3313 



G 



^1101 

^4n 1 1 



G 



1101 

42.1101 



G 



5ti2011 
g-' 22, 2011 



63 V O-t 22,3314 



1 7^0011 _ 2 /^rpOOU _ 1 / 15 7^0011 
S-' 22,4011 2T"^ ''-'22,4211 fYT 22,4213 



2 /7rpOail 
ggV 0-t 22,4413 



_ / 5 7^2011 
'*' 2T 22,2213 



2./F7n2011 

n V J-i99 991 



4 712211 
1^^^22,2213 



17712211 
90 22,2211 



1 712211 
^^-^ 22,2212 



3 -'22,2212 



2 / 2712213 
T5 V 7^22,2213 



14712212 
T5-' 22,2213' 



1112 
40.1112 



^1101 
'-^22,3303 

^1101 
"-^51,0011 

^2202 
•-^31,1112 



1101 
33.2213 



_ n /T^rplllO _ r) /T^rpllll n 2 rp\112 

— ^vZ.t22,3312 ^VZ.t22,3312 ^y 7 22,3312 

712211 _ 3^/77712211 
-'22,2211 5^1^-^22,2212 



T-.1111 I 1 T-I1112 

-^22,3313 '^ ^-^22,3313 



1 T-illl2 

^42-' 22,3314 

I 3 712212 I 
+ 5-' 22,2213 ^ 

1 T-,1110 
^-t 40,1110 

71IIOI _ t-tOOOO 
-'40,1101 -'■40,2000' 

71IIOI _ T-iOOOO 
-'42,1101 -'42,2202' 

_ 1 [zbrp^OOO _ 1 
2VT-^22,2202 2 

71IIOI I rjiOOOO 

-'51,0011 + -'51,1111' 



3_ 712212 
yjO 22,2212 



2712211 
5 -'22,2213 



71OOII I 1 / 3 712213 

6-' 22,4413 '5V 5 22,2213' 



'3t-i1111 
5-'40,llll 



^1112 
-'40,1112 



3 71OOII 
^-'40,2011' 



5t-i2202 
3-* 22, 2202 



2I71IIOI 



Tr 



22,3101 



71IIOI 
-'22,3303 



' 15t-'0000 
T"-* 22,4202' 



— _ /5t-i2000 I 1 T-12202 , ^-12202 _ 

^ Y 3-'31, 1111 I "^-'31,1111 ^" -'31,1112 



571IIOI 
3 -'31,2011 



1101 



V3^31,2211 



3t-'0000 
5 -'31,3111 



I 71IIOI 
"T -'31. 2212 



3713101 
5 -'31,0011' 



713303 I 7^/7-12202 _ /7t-i2202 _ 9 /KrpllQl , 711101 I ^0000 

-'33,0011 I V 0-133^1111 VO-t33^iii2 ^ V^-t 33 2212 + -'33,2213 ' -'33,3313' 



(C171) 



(C172) 

(C173) 
(C174) 
(C175) 

(C176) 

(C177) 

(C178) 
(C179) 



27 



^3101 

"-^11,2212 



^0011 



G' 



0011 
42,2212 



^0011 

"-r/14 99 



g; 



2011 
20,2011 



X-.2011 
Lron 99 



^2011 
Lt99 99 



5 Tn2202 
T-^11,3312 



10/712202 
T-^11,3313 



5^2000 I 1 ^2202 , ^2202 _ 

i-^11,3111 I "^-^11,3111 + -'11,3112 

1 /TTrpOOOO _ 1 [brpAOOO , 5_ti4202 , 5 t-i4202 _ 1 [b rplWl 

14 V J^o-t 11,5111 2\/3 ii'im tVI^^'^^^^ "^^^'^^^^ ^ys 11'4011 



7v^ 



7^11( 

^11, 



I rpSlOl 

"T-' 11,2212 

3 T-illlO 

i^^-t 40,1110 

^V iT ^2, 1112 + 



St^IIOI 
f-* 11,4212 " 

5 Tisaoa 

2T-'ll,2212 

S/Tnllll 

■5-'-40,llll 

^j.0011 



1 CT-SlOl 

J^o-'ii.ooii 



10rTi3303 
51^11,2213' 

1 rpOOll 
2V5^40,2011 



5 7^3101 
3 -'11,2011 



1 7^31 
^^-'ll 



3101 



/3 



42,2011 



37^0011 
5-^42,2211 



-r; 






iOOll 



42,2212 



7ti0011 
5-^42,2213' 



TnOOll _ T-T1112 
-^44,2213 -'■44,1112' 



-12211 



0^2011 , 2rTiZ 

-* 20,2011 + s-' 20,2211 



TJf-' 20,2212 



ArplllO 
■5-^20,3110 

rTi2213 

3-^20,2213' 



4^^.^1111 



-.1112 



_ 3t^0011 

20,3111 "yl -^ 20, 3112 T" 5-^20,4011 



Tn2011 2 T-rlllO 

-* 20,2211 "y5-'20,3110 



1 rp2211 
'2V5^20,2211 



6t-i1111 
5-* 22,3111 



_\_rp2212 

5^^-^20,2212 



3t'1111 
5-^20,3111 



1^-11112 
■5-'20,3112 



3 ^1112 
35 -* 20,3312 



3ti0011 
7-^20,4211 



/105" 



20.2213' 



(C180) 
(C181) 

(C182) 
(C183) 

(C184) 

(C185) 



4 T-illl2 
731 -f 22,3313 



3t'1112 
7-'22,3112 



4 rpLLLZ 

Tff-* 22,3314 



6 T-rlllO 
735^22,3312 



6 



-illll 



Tn V OJ9. 



/35-^ 22,3312 

/SrrOOll 



2 T-illl2 
TTf-t 22,3312 



5-^22,4211 



7-* 22,4212 



-.Tl 



0011 



7735^22,4213 



5t^0011 
3 -t 22,4413 



+-'22,2212 + ^7^-'^ 



2/2 rTi2212 
"5 V 35 -'22,2213 



/35" 



fT. 



1 rp2011 
273^22,2011 

2011 



67^2213 
7-' 22, 2213- 



1 rTn2011 

715^' 



22,2213 



l./^rp2212 

5 V '-'22,2212 



22,2211 " 
4 



43 Tn2211 
2573^22,2211 



25^7 



T. 



22.2213 



(C186) 
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